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Abstract 

Quantum eigenstates undergoing cyclic changes acquire a phase 
factor of geometric origin. This phase, known as the Berry phase, or the 
geometric phase, has found applications in a wide range of disciplines 
throughout physics, including atomic and molecular physics, condensed 
matter physics, optics, and classical dynamics. In this article, the basic 
theory of the geometric phase is presented along with a number of 
representative applications. 

The article begins with an account of the geometric phase for cyclic 
adiabatic evolutions. An elementary derivation is given along with 
a worked example for two-state systems. The implications of time- 
reversal are explained, as is the fundamental connection between the 
geometric phase and energy level degeneracies. We also discuss meth- 
ods of experimental observation. A brief account is given of geometric 
magnetism; this is a Lorenz-like force of geometric origin which appears 
in the dynamics of slow systems coupled to fast ones. 

A number of theoretical developments of the geometric phase are 
presented. These include an informal discussion of fibre bundles, and 
generalizations of the geometric phase to degenerate eigenstates (the 
nonabelian case) and to nonadiabatic evolution. There follows an ac- 
count of applications. Manifestations in classical physics include the 
Hannay angle and kinematic geometric phases. Applications in optics 
concern polarization dynamics, including the theory and observation 
of Pancharatnam's phase. Applications in molecular physics include 
the molecular Aharonov-Bohm effect and nuclear magnetic resonance 
studies. In condensed matter physics, we discuss the role of the geo- 
metric phase in the theory of the quantum Hall effect. 
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1 Introduction 



In a paper published in 1984, Sir Michael Berry (1984) found that a quantum 
system undergoing adiabatic evolution acquires a phase factor of purely ge- 
ometrical origin. This discovery, now called the Berry phase, or the geomet- 
ric phase, initiated a tremendous surge of research in a range of disciplines 
across physics. Berry and subsequent workers quickly established connec- 
tions to a number of apparently disparate phenomena, including amongst 
others the Aharonov-Bohm effect (Aharonov and Bohm 1959), its molec- 
ular physics analogue (Mead and Truhlar 1979), polarization optics (Pan- 
charatnam 1956), the quantum Hall effect (Thouless et al. 1982), and the 
Foucault pendulum (Eco 1989). Beyond providing a unified description of 
these phenomena, the geometric phase led to predictions of new effects, sub- 
sequently observed experimentally, and a number of substantial extensions 
and generalizations. Among these, the discovery by JH Hannay (1985) of an 
analogous effect in classical mechanics, the Hannay angle, stimulated much 
new research in purely classical physics. 

The mathematical phenomenon which underlies the geometric phase is 
holonomy. Holonomy describes transformations in a quantity induced by 
cyclic changes in the variables which control it. An example from geometry 
is the rotation of vectors parallel transported along closed paths on a curved 
surface. The mathematics of holonomy and curvature has long played a 
role in physics through general relativity and quantum field theory. With 
Berry's work came the appreciation of a broader scope of application. When- 
ever a physical system is divided into parts and one attempts to describe a 
subsystem in isolation, the influence of the others is manifested through ge- 
ometric phase effects. The pervasiveness througout physics of this reductive 
procedure accounts for the ubiquity of the geometric phase. Introducing 
a collection of major papers in the field, A Shapere and F Wilczek(1989) 
wrote, "We believe that the concept of a geometric phase, repeating the 
history of the group concept, will eventually find so may realizations and 
applications in physics that it will repay study for its own sake, and become 
part of the lingua franca." 

An introductory account of the geometric phase is given in Section 2. 
This is followed in Section 3 by a discussion of several theoretical develope- 
ments. Applications to classical systems, optics, and molecular and con- 
densed matter physics are given in Sections 4-6. The topics chosen for 
discussion meant to be representative rather than comprehensive. Among 
several omissions, we mention the theory of fractional statistics, which is 
covered extensively in the volume by Shapere and Wilczek (1989a), and ap- 
plications to elementary particle physics, which are reviewed by Aitchison 
(1987). The historical development of the subject is discussed by Berry 
(1990a). 
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2 Basic Account of the Geometric Phase 



The geometric phase was discovered by Sir Michael Berry in 1983 in the 
course of a critical reexamination of the adiabatic theorem in quantum me- 
chanics. This theorem says, in essence, that if a quantum system is changed 
sufficiently slowly, it does not undergo transitions. Thus a system initially 
in an eigenstate of a slowly changing Hamiltonian remains in the evolved 
eigenstate. Berry considered the overall phase accumulated by the changing 
eigenstate and found an unexpected contribution. This extra phase is of 
a geometrical nature, in that it does not depend on the rate at which the 
adiabatic process is carried out. 

One way that a slowly varying Hamiltonian can arise physically is when 
the system under consideration interacts with another, much slower system. 
This leads to the consideration of the quite general problem of two coupled 
systems evolving on different time scales. In this context, the geometric 
phase describes an effect on the fast system produced by the slow one. As 
action in physics is accompanied by reaction, there is a reciprocal effect on 
the slow system generated by the fast. This takes the form of an effective 
vector potential in the slow Hamiltonian; the associated Lorentz-like force is 
called geometric magnetism. The same geometry underlies both geometric 
magnetism and the geometric phase. 

This geometrical structure is intimately connected to degeneracies in 
the energy spectrum. A typical Hamiltonian has a nondegenerate spectrum 
(excepting systematic degeneracies due to symmetry, which are discussed 
in Section 3.2). Varying a single parameter (eg time itself in the case of a 
time-dependent Hamiltonian) is unlikely to produce any degeneracies; typi- 
cally two or three parameters (the number depends on whether time-reversal 
symmetry is present or not) must be adjusted independently in order to find 
one. Nevertheless, the presence of degeneracies in Hamiltonians which are, 
in some sense, near the Hamiltonians of the adiabatic process, are largely 
and sometimes wholly responsible for the geometric phase. This connection 
underlies several applications, including to the quantum Hall effect (Sec- 
tion 6.3). 

2.1 Quantum adiabatic theorem 

Adiabatic invariants are quantities which are nearly conserved when the 
environment of an isolated system is slowly changed (slowly means in com- 
parison with the system's internal dynamics). Adiabatic theorems provide 
estimates as to how accurately this conservation is respected. In quantum 
mechanics, the adiabatic invariants are the probabilities P n = |(n|^)| 2 to 
be in energy eigenstates. These are constant for stationary Hamiltonians, 
but vary if H is time-dependent. The quantum adiabatic theorem states 
that as dH/dt goes to zero, so do the probabilities of transitions between 
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eigenstates (the rate depends on how smooth dH/dt is), so that the P n are 
nearly constant. Therefore, within the adiabatic approximation, dynamics 
under a slowly varying Hamiltonian H{t) is essentially no more complicated 
than under a stationary one. In particular, if \ip) begins in an eigenstate, it 
remains in one ( nonadiabatic transitions are discussed in Section 3.5), and 
its evolution is determined up to a phase factor. 

We shall now examine this phase. Suppose that at t = 0, \tp) is in the nth 
eigenstate |n) (assumed nondegenerate) of H. For its subsequent evolution 
we may write 

|Y>(t)) » exp f-i/h J E n {r)dT^j \n{t)). (2.1) 

In so doing we have explicitly identified part of the overall phase, the time 
integral of the frequency E n /h, which is called the dynamical phase. It is 
just the phase one would expect for a solution of the Schrodinger equation. 
Any remaining phase is contained implicitly in the eigenstate \n(t)). Its 
phase is determined by requiring that (n\(ihdt — H)\tp) vanish (it would 
automatically if \ip(t)) were an exact solution of the Schrodinger equation). 
This leads to the condition 

(n\n) = 0, (2.2) 
which determines the phase of \n(t)) given an initial choice at t = 0. 



2.2 The geometric phase 

What Berry discovered is that condition (2.2) yields a path-dependent phase 
factor. This phase factor is particularly interesting when the Hamiltonian 
is cycled, ie made to return to itself after some time. Then it depends only 
on the cycle itself, and not, in contrast to the dynamical phase, on the rate 
at which the Hamiltonian is cycled. For this reason it is called, as well as 
the Berry phase, the geometric phase. 

To discuss this further, it helps to formulate the problem slightly differ- 
ently. Suppose the time dependence of the Hamiltonian is expressed through 
its dependence on some external parameters R = (X, Y, Z), which are made 
to vary slowly in time along a path R(i). Then the eigenstates |n(R)) de- 
pend on R as well, and (assuming no degeneracies) are determined for each 
R up to an overall phase. Let us choose this phase to vary smoothly but 
otherwise arbitrarily with R. By fixing the phase independently of the dy- 
namics, we may no longer satisfy the condition (2.2) along every path R(t). 
We therefore allow for an additional phase factor exp(ry n (t)) in the adiabatic 
ansatz (2.1), which becomes 

\i/>(t)) « exp E n (R{r))dT^ exp(i 7n (t))|n(R(*))). (2.3) 
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The condition {n\(ihdt — H)\ip) = now yields an evolution equation for j n . 
This may be written in the form j n (t) = — A n (R(t)) • R(i), where A n (R) is 
a vector field in parameter space given by 

A„(R) = Im(n(R)|Vn(R)>, (2-4) 

or equivalently, in terms of eigenfunctions ip n (x, R) = (x\n(R)}, by A„(R) = 
Im j ip n (x, R)V , 0(x, R)cfe (here V acts on the parameters R). 

If the Hamiltonian is cycled around a closed curve C in parameter space, 
then \ip(t)) acquires a geometric phase 

7„(C) = — [ A n -Kdt = <f A n ■ dR. (2.5) 
Jo Jc 

Its expression as a line integral reveals its geometrical character; the geo- 
metric phase depends only on the cycle C, and not on the manner in which 
the cycle is traversed in time. Using Stokes' theorem, 7„ can be expressed 
as 

7n(C) = - / V n (R)-dS, (2.6) 
Js 

the flux of the vector field V„ = V x A n through a surface S in parameter 
space whose boundary is C. V n is given by 

V n = Im(Vn| x |Vn). (2.7) 

From (2.6) it is evident that only cycles which enclose nonzero area can 
have nontrivial geometric phases. For example, j n (C) vanishes for cycles 
consisting of a path and its retrace. 

It is easy to see that the geometric phase does not depend on the choice 
of phase for |n(R)). Indeed, under the change of phase 

KR))->exp(i x (R))|n(R)>, (2.8) 

A n (R) transforms according to A n (R) — > A„(R) + Vx(R); this shift by a 
perfect gradient leaves its integral round a closed loop, 7„, and its curl, V n , 
unchanged. Thus, the geometric phase is an intrinsic property of the family 
of eigenstates |n(R)). 

It is useful to draw an analogy with electromagnetism, in which A n 
corresponds to the vector potential and V n to the magnetic field. The geo- 
metric phase 7„ corresponds to a magnetic flux. The change of phase (2.8) 
corresponds to a gauge transformation, which shifts the vector potential by 
a gradient and leaves the magnetic field unchanged. 

Substituting the first-order perturbation expansion 

|Vn) = £ <MVH\n) + (n|Vn) |n) (2 g) 



8 



into (2.7) leads to a useful formula, 



v-^ (n\VH\m) x (m\VH\n) 

In this form, the independence of V n on the eigenstate phases is manifest. 
Another formula with this property is 

V n = Tr {p n Vp n X Vpn), (2.11) 

where p n = \n)(n\ is the density operator. It too is useful in applications. 

One of the applications discussed by Berry (1984) is to the Aharonov- 
Bohm effect (Aharonov and Bohm 1959). The Aharonov-Bohm effect de- 
scribes the dependence of quantum mechanics on electromagnetic fields in 
regions of space which are physically excluded. This dependence is medi- 
ated through the vector and scalar potential, whose relation to the electric 
and magnetic fields is nonlocal. The most celebrated example concerns an 
electron moving in the vicinity of an inpenetrable perfect solenoid of infinite 
extent. Since outside the solenoid the magnetic field is zero, the classical me- 
chanics of the electron is insensitive to any magnetic flux inside the solenoid. 
But the quantum mechanical behaviour of the electron does depend on the 
interior flux, as Y Aharonov and D Bohm showed by calculating its scat- 
tering cross-section (a prediction subsequently experimentally confirmed by 
Chambers (I960)). 

The Aharonov-Bohm effect can be understood in terms of interference 
between de Broglie waves passing on opposite sides of the solenoid. The 
phase difference between the two waves includes the integral 

7AB = Y c f c A(s) ( R )- dK ( 2 - 12 ) 

of the vector potential A( S )(R) of the solenoid round a loop C enclosing it. 
(In (2.12), c is the speed of light and — e the electron charge.) The Aharonov- 
Bohm phase ^ab is just 2tt times the magnetic flux in the solenoid measured 
in units of the flux quantum hc/e. 

It may be regarded as a geometric phase in the following manner. We 
imagine the electron confined to a large box, centred at R, which is not 
penetrated by the solenoid. The eigenstates of the electron |n(R)) depend 
parametrically on R; changes in R induce gauge transformations in the 
eigenfunctions. One finds that Im (n(R)| Vn(R)) = e/hcA^ s \ so that in 
this case the geometric vector potential (2.4) coincides with the physical 
vector potential of the solenoid. On being transported (not necessarily adi- 
abatically in this case) round the solenoid, the electron eigenstates acquire 
a geometric phase which is precisely the Aharonov-Bohm phase. 
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2.3 Two-state systems 

Two-state systems provide a canonical example of the geometric phase. The 
requisite calculations are easily carried out, and the main features of the 
general case are represented. The parameterized family of Hamiltonians can 
be written in the form 

H(R) = E(K)I + F(R) • <r, (2.13) 

where I is the identity and a are the Pauli spin matrices. The term E(R)I 
shifts the energy by a constant and affects neither the eigenstates nor their 
geometric phases. Both of these are determined by the field F(R). 

An illustrative example is a spin- 1/2 particle in a uniform magnetic 
field. The parameters R are then just the components of the field, so that 
F(R) = ^R, where fi is the magnetic moment. The eigenstates 

i+f^U^/O- i-< r »K«0 (2 - i4) 

depend on the polar angles (0,<fi) of R, and the energy levels £±(R) = 
±^/j,hR on its magnitude. Straightforward calculation using (2.4) and (2.7) 
yields 

Thus, A± and V± are respectively the vector potential and magnetic field of 
a magnetic monopole of strength ±2ir at the origin. Under the influence of 
a slowly varying magnetic field R(i) taken through a cycle C, a spin initially 
in an up-state returns to itself with a dynamical phase —fi/2h f Rdr and a 
geometric phase 7+(C) = — f s V + • dS. Gauss's law implies that 

7+(C) = — i x solid angle subtended by C. (2-16) 

This spin- 1/2 geometric phase has been observed in neutron spin rotation 
experiments by Bitter and Dubbers (1987). A special case is when C is 
confined to a plane containing the origin. Then its solid angle is 2ir times 
the number of times m the origin is enclosed, and the corresponding phase 
factor is (— l) m . This corresponds to the well-known sign change of spin-1/2 
particles under rotations through 2tt. 

It is straightforward to generalize these results to the case of arbitrary 
F(R) in (2.13). The flux field is then given by 

V ± (R)^4„ /' V ^ VFt (2-17) 

(eijk is the antisymmetric Levi-Cevita symbol), and it has monopole sources 
at points where F vanishes. The monopoles have strength ±2ira, where a is 
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the sign of the (assumed nonvanishing) determinant det dFi /dRj evaluated 
where F = 0. 

As is familiar from Dirac's analysis of magnetic monopoles (Dirac 1931), 
there is necessarily a string of singularities of the vector potential attached 
to the monopole. The position of the string depends on the choice of gauge, 
and Dirac's quantization condition for electric charge follows from demand- 
ing that these choices be physically indistinguishable. In the context of the 
two-state geometric phase, the monopole strings of A±(R) lie along disconti- 
nuities in the phases of |±(R)). When the eigenstates are chosen as in (2.14), 
these discontinuities lie along the negative z axis. By varying the phase of 
the eigenstates (as in (2.8)), the string can be made to pass in and out of 
a given cycle C. As it does, the geometric phase 7±(C) changes by ±2ir. 
Therefore, the geometric phase factor exp(z7±(C)) remains unchanged, and 
there are no observable consequences. 

2.4 Degeneracies 

In the spin-1/2 example, the flux fields V±(R) are generated by point 
sources of strength ±2ir at the origin R = 0. What makes the origin spe- 
cial is that it is the only point in parameter space where the energy levels 
are degenerate. This relationship between the geometric phase and energy 
level degeneracies turns out to be quite general. In the absence of symme- 
tries, a typical Hamiltonian in a parameterized family has a nondegenerate 
spectrum. However, at certain points in parameter space, belonging to the 
so-called degenerate set, two (or more) energy levels become degenerate. In 
the absence of time-reversal symmetry (the time-reversal symmetric case is 
discussed in the following section), the dimension of the degenerate set is 
typically three less than the dimension of parameter space itself. Thus, for 
the three-parameter families H(R) we have been considering, degeneracies 
occur at isolated points. We have already seen this to be the case in the 
spin-1/2 example. (Note that there, the magnetic field is responsible for 
breaking time-reversal symmetry.) 

More generally, let us consider the behaviour of the flux field V n (R) in the 
neighbourhood of a degeneracy R* between two eigenstates, say the nth and 
n— 1st, of an arbitrary family of Hamiltonians. As E n (R) — E n -i(R) is small 
near R*, perturbation theory (cf (2.9)) implies that the dominant compo- 
nent of | Vn) lies along \n— 1). Therefore, for the purposes of approximating 
V n (which is defined in terms of |Vn)), the Hamiltonian can be restricted 

to the nearly-degenerate subspace spanned by the states |1) d = n |n(R*)) and 

1 2) d = n \n — 1(R*)). One obtains in this way a two-state system of the form 
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(2.13), where the field F(R) is given by 



F x = Re (l|i?(R)|2), F y = lm(l\H(R)\2), F z 



(l\H(K)\l) - <2|tf(R)|2) 
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(2.18) 



In the neighbourhood of R*, V ra (R) is given approximately by (2.17). 
There it has a single monopole source at the degeneracy R* of strength ±2ir. 
We may apply a similar analysis in the neighbourhood of each degeneracy 
of \n) to obtain 



Thus, the flux field for the nth eigenstate has monopole sources of strength 
±2ir at points where \n) is degenerate. (Their sign depends on whether the 
degeneracy is with the state above or below, and also on the sign of the deter- 
minant det dFi (R*)/dRj.) V n may in addition have a purely divergenceless 
contribution. 

It follows that the flux of V n through a closed surface in parameter space 
is 2ir times the number of degeneracies, counted with their signs, contained 
within. 

2.5 Time-reversal symmetry 

The presence of time-reversal symmetry imposes constraints on the geomet- 
ric phase. These are strongest when the eigenstates themselves are time- 
reversal symmetric, so we consider such cases first. Then the wave functions 
Y> n (x, R) = {x\n(R)) can be chosen to be real, and the geometric phase 
factor is necessarily a sign ±1, so that 7„ is an integer multiple of n. In this 
case, the geometric phase depends only on the topology of the closed curve, 
not on its geometry. As C is deformed, unless J n (C) changes discontinu- 
ously, it necessarily remains constant. It can change (discontinously) when 
C passes through degeneracies (points R* where |n) is degenerate). For sys- 
tems with time-reversal symmetry, the dimension of this degenerate set is 
typically two less than the dimension of parameter space, so that it consists 
of lines in three-dimensional parameter space, or points in two-dimensional 
parameter space. The geometric phase j n {C) is ir times the number of de- 
generacies enclosed by C, counted with the appropriate sign. An important 
consequence is that the flux field V n has delta-function singularities at the 
energy level degeneracies and is zero elsewhere. 

The spin-1/2 example of Section (2.3) provides an illustration if we fix Y 
(the y-component of the magnetic field) to be zero. Then the spin Hamilto- 
nian i^B.tr is real symmetric rather than complex Hermitian, and we can 
take its eigenvectors to be real. The parameter space is the XZ plane, and 
the solid angles subtended by curves in this plane are multiples of 2ir. The 
geometric phase 7±(C) is just ±ir times the number of times C encircles the 




(2.19) 
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origin. Another example concerns the electronic ground state of triatomic 
molecules A3 (cf Section 6.1). 

In cases where the eigenstates are themselves not time-reversal symmet- 
ric, and instead occur in pairs related by time-reversal, the constraints on 
the geometric phases are much weaker. They are no longer necessarily mul- 
tiples of 7r (ie, the phase factors are no longer necessarily real). However, 
geometric phases for states related by time-reversal necessarily have opposite 
signs. 

2.6 Experimental observation 

Since the physical properties of a quantum system does not depend on its 
overall phase, one could wonder whether the geometric phase is observable. 
Numerous experiments have shown that it is. 

Most direct experimental studies of the geometric phase belong to one 
of three types. The first is the 'one state, two Hamiltonian' experiment. 
An ensemble of states initially in an energy eigentate is coherently divided 
into two streams. The first evolves under a constant Hamiltonian, while the 
Hamiltonian driving the second is cyclically varied. The two streams are 
then recombined and allowed to interfere, and their relative phase deter- 
mined. The relative phase includes both dynamical and geometric contribu- 
tions, and to extract the geometric term the dynamical phase must either be 
made to vanish, or else subtracted off explicitly. Examples are the optical 
polarization experiments discussed in Section 5.2. 

The second type is the 'two state, one Hamiltonian' experiment. Here 
an initial state is prepared in a superposition oti\m) + fti\n) of energy 
eigenstates. The Hamiltonian governing its evolution is taken through an 
adiabatic cycle. The final state |Y>/) = cn/|m) + Pf\n) is physically dis- 
tinguishable from the initial one, as the relative phase of the superposed 
eigenstates has changed. Part of this (observable) change in the relative 
phase involves the difference 7 m — 7„ in geometric phases. The observa- 
tion of optical activity in coiled fibres by Tomita and Chiao (1986) (see 
Section 5.1) provides an example of this second type. 

In the third type of experiment, the adiabatic cycle is performed not 
once but repeatedly, with period T. With each cycle, |n) acquires a phase 
4> = — f E n (R,(t)dt/h + 7„, so that there appears in its mean angular fre- 
quency oj = (f>/T a geometric shift 5oj g = 7 n /T. This frequency shift can be 
detected in spectroscopic measurements, and is the basis for nuclear mag- 
netic resonance studies of the geometric phase (cf Section 6.2). 

2.7 Geometric magnetism 

So far we have been regarding the evolution of the parameters R as exter- 
nally prescribed. Now let us suppose the parameters of the Hamiltonian 
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H(R) are dynamical variables in their own right. For example, they could 
be the coordinates of a heavier (and therefore slower) system to which the 
original one, the fast system, is coupled. For definiteness, we take the Hamil- 
tonian of the coupled system to be of the form P 2 /2M + H(R), where P 
is the momentum of the slow system, and the dependence on the fast co- 
ordinates and momenta is left implicit in H(R). Both systems are to be 
treated quantum mechanically. The geometric phase describes an effect on 
the fast system due to the slow. There is a corresponding influence on the 
slow system due to the fast. 

This is the natural setting for the Born-Oppenheimer approximation, in 
which the fast and slow systems are the electrons and nuclei of a molecule. 
(For consistency with the preceding discussion, we denote the nuclear co- 
ordinates by a single vector R; it is easy to adapt the notation to accomo- 
date more than one nucleus.) By assuming that the electronic state rapidly 
adjusts to changes in the nuclear configuration, one is led to look for ap- 
proximate eigenfunctions of the form * n (r, R) = ^ ra (R)^ n (r, R), where the 
electronic wavefunction (j) n (r, R) is an eigenstate with energy E n (H) of the 
electronic Hamiltonian H(R), in which the nuclear coordinates appear as 
parameters. The nuclear wavefunction ifr n (R) is an eigenstate of the nuclear 
Hamiltonian — /i 2 V 2 /2M + E n (R); the electronic energy plays the role of 
a potential. The Born-Oppenheimer approximation consists of neglecting 
the weak variation of the electronic wave function (f) n (r, R) with the nuclear 
coordinates. (It is essentially the adiabatic approximation (2.1) adapted to 
coupled systems.) 

This was the context for seminal work in this subject by CA Mead and 
D Truhlar in 1979. The primary objective of this work, which anticipated 
certain aspects of the geometric phase, was to derive corrections to the Born- 
Oppenheimer approximation. Such corrections appear when the momentum 
operator P = — i/iV is applied to the putative eigenstate ^ n (R)|n(R)). The 
result, —ihCVifj n \n) + ip n \ Vn)), is no longer proportional to the electronic 
eigenstate |n(R)). However, if we neglect electronic components orthogonal 
to |n), we obtain — zft(Vy> n |n) + ?/> n (n|Vn}) |n), which is. The effect on the 
nuclear Schrodinger equation is to redefine the momentum operator, 

d i n (-ihV + ftA n )Y>, (2.20) 

Here A„(R) is given by (2.4). Thus, the field whose line integral gives 
the electronic geometric phase appears as a vector potential in the nuclear 
Hamiltonian. There it corresponds to an effective magnectic field V„ = 
V x A n which modifies the nuclear dynamics. This is called geometric 
magnetism. 

The situation described above is quite general. Whenever there are two 
coupled systems evolving on different time scales, we can effect an approx- 
imate description of the slow system by averaging over the fast motion. 
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There appears in the slow dynamics an effective vector potential and asso- 
ciated magnetic field. The vector potential and magnetic field are precisely 
the fields whose line integrals and fluxes describe the geometric phase. (In 
cases where the coupling involves not only the slow coordinates but the slow 
momenta as well, this description has to be modified somewhat.) 

Geometric magnetism is responsible for a variety of phenomena. Mead 
and Truhlar analysed its effects in nuclear rotation- vibration spectra, which 
have since been observed experimentally (cf Section 6.1). These effects can 
be nontrivial even when time-reversal symmetry is present (cf Section 2.5), 
due to the change in sign of the electronic wavefunction (f> n (r, R) when the 
nuclear coordinates R are carried round a degenerate configuration. Like 
the magnetic field of the Aharonov-Bohm effect (Section 2.2), the flux field 
V„(R) in this case vanishes everywhere except at these degenerate config- 
urations, where it is singular. Mead (1980) has called this the "molecular 
Aharonov-Bohm effect" . 

Geometric magnetism survives in the classical limit. When this limit is 
applied to the slow system alone, it leads to geometric phase modifications 
of the semiclassical (Bohr-Sommerfield) quantization conditions. These are 
discussed by Kuratsuji and Iida (1985), Wilkinson (1984), and Littlejohn 
and Flynn (1991). 

The situation where both the fast and slow systems are fully classical 
is discussed in Section 4.1.5. Here let us mention a simple but compelling 
illustration discussed by Li and Mead (1992). Let us consider the Born- 
Oppenheimer approximation applied to an atom, rather than the molecule. 
The electronic energy levels E n (R) depend on the atom's centre-of-mass R, 
and appear as potentials in the equations of motion R = — V-E n + F(R), 
along with the external force F(R). In the absence of external potentials, the 
energy levels are independent of R, and F = 0; the centre-of-mass behaves 
as a free particle. Now suppose the atom moves in a uniform magnetic field. 
The energy levels are still R-independent (by translation invariance), but 
there now appears in the centre-of-mass dynamics the Lorentz force F = 
+ZeR x B. The predictions of this analysis are spectacularly unphysical; 
walking through the earth's magnetic field, you would experience a vertical 
force sufficient to hurl you into the air. What is missing from this analysis 
is geometric magnetism from the electrons; the field V n (R) introduces into 
the centre-of-mass dynamics a second Lorentz force which either exactly 
or nearly (in case the electron eigenstate has a nonzero magnetic moment) 
cancels the first. 

3 Theoretical Developments 

A number of theoretical developments of the geometrical phase are de- 
scribed. These include its mathematical setting in terms of holonomy under 
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parallel transport, as well as several generalizations which remove the as- 
sumptions of nondegeneracy and adiabaticity. 

3.1 Geometric phase as holonomy 

Objects rigidly transported round a closed path in a curved space can return 
with a different orientation - this is holonomy. Here is a simple (and surely 
physical) example. Hold your right arm outstretched before you with your 
palm facing down. Keeping your wrist rigid, swing your arm into your chest, 
then raise it clockwise in front of you, and finally swing it away downwards 
until it is outstretched before you again. Your arm has returned to its 
original position, but your hand has turned clockwise through 90°. 

The mathematical setting for this phenomenon is the theory of connec- 
tions on fibre bundles. What follows in Sections 3.1.1— 3.1.4 is an informal 
account of the relation between this theory and the geometric phase. No at- 
tempt is made at mathematical precision; the intention is to convey a sense 
of the principal ideas involved. Some standard references are (Eguchi et al 
1980) and (Nash and Sen 1983). 

3.1.1 Fibre bundles 

A fibre bundle E is a space, or manifold, composed of two smaller ones, 
the base manifold B and the fibre F. The fibre bundle is constructed by 
attaching a copy of the fibre to each point of the base in a particular way. 
The simplest trivial fibre bundle, is where E is just the cartesian 

product B x F of the base with the fibre. A general fibre bundle cannot be 
so expressed. While it is always possible to partition it into smaller pieces 
which are themselves cartesian products M a x F of subsets of the base with 
the fibre, these pieces may fit together in an interesting and nontrivial way, 
which reflects the topological properties of the bundle. 

The configuration space E of an outstretched arm can be described as a 
fibre bundle. In this case, the base manifold is the unit sphere, whose polar 
coordinates (9, 4>) determine the arm's direction. The fibre is the unit circle, 
whose single angle coordinate x describes the orientation of the hand relative 
to some reference position. In this case E is not simply the cartesian product 
M x F. This is because it is not possible to assign, in a continuous way, a 
reference hand position for every arm direction. This is a consequence of the 
"hairy ball theorem", ie the fact that it is impossible to comb flat a hairy 
ball without introducing singularities (eg crowns, cowlicks, bald patches). 

As in the arm example, it is often the case in physical applications that 
the fibre describes internal degrees of freedom. Usually the fibre has some 
additional mathematical structure as well. If it is a real or complex vector 
space (describing spin degrees of freedom, for example), then E is called a 
real or complex vector bundle. If the vector space is one-dimensional, E is 
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called a line bundle (real or complex). Or else, the fibre might be a group 
(for example, the two- or three-dimensional rotation group, describing the 
axial or spatial orientation of a rigid body). Then E is called a principal 
bundle. 

3.1.2 Connections and curvature 

To a given path P on the base there correspond innumerable paths in the 
fibre bundle which project down to it; a motion in the base leaves the fibre 
motion undetermined. A connection is a rule for associating to P a particular 
path P through a given point in the bundle, the lifted path (see Fig. 1). A 
point moving along the lifted path P is said to be parallel transported. The 
connection rule is expressed in differential form by making the fibre velocity 
a linear function of the base velocity, which can also depend on the base and 
the fibre positions. For vector bundles, the dependence on fibre position 
should be linear as well. There are analogous restrictions in the case of 
principal bundles. In physical applications, a particular connection is often 
suggested by kinematical or dynamical considerations, eg, moving your arm 
without unnecessarily twisting your wrist. 

Let us consider the lift C of a closed path C on the base (see Fig. 2). 
Its endpoints necessarily lie in the same fibre, but they need not coincide. 
If they do not, the difference between the endpoints is called the holonomy 
of C. In the outstretched arm example, C is the spherical right triangle 
described by your arm, and its holonomy is the ir/2 turn of your hand 
which results from traversing it in the manner prescribed. In case there is 
holonomy for arbitrarily small closed paths, the connection is said to have 
curvature. Curvature describes the infinitesimal holonomy acquired round 
an infinitesimal closed paths, and may be expressed in terms of derivatives 
(a generalized curl) of the connection. 

Let us consider a very simple example. We take the base M to be three- 
dimensional Euclidean space {R G M 3 }, the fibre F to be one-dimensional 
complex vector space {z G C}, and the bundle E to be their cartesian 
product {(R,z) G R 3 x C}. Thus E is a trivial complex line bundle. A 
connection specifies the fibre velocity z as a function of the fibre position 
z, R and R whose dependence on z and R is linear. This relation may 
be expressed as z = (W(R) ■ R) z, where W(R) is a complex vector field. 
W(R) is called the connection form, or simply the connection. 

Given a motion R(i) in the base, its lifted motion in the fibre is given by 

z{t) = exp ^ W(R) • dRj zq, where zo denotes the initial position in the 

fibre of the lifted path. If R(t) describes a closed path, then its holonomy is 
the scalar factor expT, where T = W(R) • dR. Using Stokes' theorem, 
we may express T as the flux of a vector field V = V x W through a surface 
S bounded by C. V is the curvature of the connection W. 
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Of particular relevance to the geometric phase is the case in which the 
connection preserves the length \z(t)\ of vectors along the lifted path. Then 
the connection is said to be hermitian. The connection form W(R) is purely 
imaginary, and may be written as zA(R), where A(R) is real. 

3.1.3 The geometric phase revisited 

After seeing a preliminary version of Berry's 1984 paper, B Simon (1983) 
observed that the geometric phase finds its natural mathematical expression 
as the holonomy of a connection on a fibre bundle. The fibre bundles in 
question are the complex line bundles associated with nondegenerate 
eigenstates of a quantum Hamiltonian H(H). The base manifold is the 
parameter space (we continue to assume this to be R 3 ), and the fibres are 
vectors in Hilbert space which satisfy the Schrodinger equation (H (R) — 
E n (R))\ip) = 0; these are just the unnormalized nth eigenstates, and are 
determined up to a complex scalar. 

Given an eigenstate \tp) of H(R), a connection determines an eigenstate 
\ip) + \dip) °f the nearby Hamiltonian H(R+dR) by fixing its norm and phase. 
On physical grounds it is sensible to require that the connection conserve 
probability, ie that \ip) + \dip) and \ip) both have unit norm. This implies 
that Re (tp\dtp) = 0. The simplest rule for fixing the phase of \tp) + \dip) is 
to take lm (ipldTp) = as well. This determines completely a connection on 
£;( n ) , which is defined by the relation 

M<ty> = 0. (3.1) 

Under parallel transport round a closed path C, an eigenstate \ip) acquires 
a phase factor exp(«7 n ). This phase factor, the holonomy of C, is precisely 
the geometric phase factor. This follows by noting that the connection (3.1) 
is just the condition (2.2) derived from adiabatic evolution. 

There is an exact mathematical analogy between the outstretched arm 
and a variant of the spin example of Section 2.3. If a spin-1 (instead of a 
spin-1/2) particle, initially in its down-state, is driven by a magnetic field 
which is slowly cycled round a closed path C, then its geometric phase 71 is 
precisely the turn in the hand of an arm describing the same circuit. Both 
are equal to the solid angle subtended by C. 

3.1.4 Chern classes 

Some fibre bundles are trivial cartesian products M x F, and others are 
not. Topology lies at the root of this distinction, and characteristic classes 
make this distinction precise. They describe an intrinsic twistedness in fibre 
bundles which cannot be combed away. Bundles whose characteristic classes 
are the same can be smoothly transformed into one another; bundles whose 
characteristic classes are different cannot. 
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For complex vector bundles, the characteristic classes are called Chern 
classes. The number of Chern classes is equal to the dimension of the fibre. 
The rth Chern class associates an integer c r (N), the Chern number, to every 
2r-dimensional, boundaryless submanifold N of the base manifold M. Thus, 
the first Chern class assigns an integer ci(S) to every closed surface S in M. 
There are various ways to calculate Chern numbers. Among these, there are 
formulas in terms of integrals over N involving the curvature. 

An illustrative example is the celebrated Gauss-Bonnet theorem from 
classical differential geometry. The theorem states that the integral f s KdS 
of the Gaussian curvature K over a closed one-sided surface S is equal to 
27r times an integer x, the Euler characteristic. The Euler characteristic is 
a topological feature of the surface. It is related by the formula x = 2 — 2g 
to the genus g, the number of handles on S when it is regarded through 
a topologist's eyes as a many-handled sphere. The Gauss-Bonnet theorem 
can be verified immediately for the unit sphere (K = 1, g = 0). What is 
remarkable is the implication that J s KdS = Air for any surface smoothly 
deformable to a sphere, but f s KdS / 4ir for any surface (such as a torus) 
which is not. 

The Gauss-Bonnet theorem can be formulated in terms of fibre bundles. 
The bundle consists of the surface S (the base) together with its tangent 
vectors at each of its points (the fibres). A connection is defined through 
the Riemannian rule for parallel transport on a curved surface. By iden- 
tifying two-dimensional real tangent vectors with one-dimensional complex 
ones, the fibre bundle may be regarded as a complex line bundle, and its 
first Chern number c\ (S) turns out to be twice the integral f s KdS of the 
Gaussian curvature. 

Returning to considerations of the geometric phase, let us consider the 
complex line bundles E^> of quantum eigenstates, discussed in the preceed- 

(n) 

ing section. Each has a single Chern class c\ which assigns an integer to 
every closed surface S in parameter space. This integer is given by 



where in the second equation we have used the divergence theorem; B is 
a three-dimensional region bounded by S. From (2.19), the last integral 
is precisely the number of energy level degeneracies (points R* £ B where 
E n (R*) = E n ±i) contained in S, counted with their appropriate sign (see 
Section 2.4). In terms of the electromagnetic analogy (Section 2.2), it is 
the number of magnetic monopoles in S counted with their charge. Thus, 
energy level degeneracies manifest themselves in the nontrivial topology of 
the eigenstate line bundles E^ N \ 




(3.2) 
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3.2 Degenerate eigenstates and nonabelian geometric phases 

If a system possesses a sufficient degree of symmetry, its energy levels will 
be degenerate. Two examples are the angular momentum multiplets of rota- 
tionally symmetric systems, and Kramers degeneracy in systems with time- 
reversal symmetry and an odd number of fermions. If the symmetry persists 
throughout the parameterized family H(R), then so do these degeneracies. 
Under a slowly varying Hamiltonian H(R(t)), the adiabatic theorem implies 
that transitions to states outside a degenerate multiplet are small. However, 
transitions to states within the multiplet need not be, no matter how slow 
the variation (transitions between degenerate states cost no energy). 

F Wilczek and A Zee (1984) generalized Berry's analysis to degenerate 
eigenstates. They showed that under adiabatic cycling round a closed circuit 
C, such states need not return to themselves up to a phase, but instead can 
evolve into new states within the multiplet. The transformation from initial 
to final states is generated by a unitary matrix U(C), called the nonabelian 
geometric phase. Nonabelian geometric phases can be interpreted as holon- 
omy, but of a more general sort than was considered in Section 3.1.3. There 
is a corresponding generalization of geometric magnetism. 

The analysis, which runs parallel to the nondegenerate case, proceeds as 
follows. Let E(R) denote an r-fold degenerate energy level whose degenerate 
eigenspace is spanned by an orthonormal basis \j(R)),j = 1,2,... ,r. We 
consider evolution under a slowly changing Hamiltonian H(R(t)) of states 
\ij}j{t)) initially equal to |j(R(0))). Neglecting transitions to states outside 
the multiplet, we look for approximate solutions of the form 



ie, a superposition of evolving degenerate eigenstates multiplied by an overall 
dynamical phase factor. For the inner products (tpk(t)\ipj{t)} to be conserved 
in time, the rxr matrix U(t) must be unitary. Its evolution is determined by 
requiring (k(R)\(ihdt— H)\if)j (t)) to vanish for each k (it would automatically 
if ipj(i) were an exact solution of the Schrodinger equation). This leads to 
the equation of motion 



where A(R) = (A^, A^, A^)(R) is a vector of r x r Hermitian matrices 
whose matrix elements are given by 




(3.3) 



U(t) 



iU(t) A(R(i)) • R = -iU(t)A(t), 



(3.4) 



A<°>(R) = -i(*(R)|V«j(R)>, 



(3.5) 



and A(t) = A(R(t)) • R(t). 
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The solution of (3.4) is given by the time-ordered product 

U(t) = T exp f-i J A(r) d^j 

= 1+1 A(r) dr + [ [ A(r)A(r') drdr 1 + ■■■ , (3.6) 

J0<T<t J J0<T<T'<t 

where I is the identity matrix. The time ordering indicates that in expanding 
the exponential of the integral, all matrices are to be ordered with later times 
to the right. If the matrices A(r) at different times commute with each other, 
then the ordering of the factors is irrelevant, and the r.h.s. of (3.6) simplifies 

to I + J 0<T<t A(r) dr + (J Q<T<t A(r) dr) 2 /2\ + • • • = exp (-i f* A(r) dr) , eg 
the ordinary exponeniated integral. In particular, if r = 1, then A(r) is a 
scalar, and U (t) is simply a phase factor. 

Suppose the Hamiltonian is taken round a closed cycle C. Then the re- 
lation between initial and final states is described, apart from the dynamical 
phase factor, by the unitary matrix 

U(C) = Texp A(R)-dR^). (3.7) 

U (C) is the nonabelian geometric phase. The terminology derives from the 
fact that the factors A(R) -dR in the time-ordered product do not in general 
commute. They do, however, in the nondegenerate case r = 1; then U(C) 
is just the ordinary (abelian) geometric phase factor exp(ij(C)). Like the 
ordinary geometric phase, U (C) depends only on the path C and not on the 
rate at which it is traversed. 

For infinitesimal circuits, U(C) is of the form I + iY(C), where T(C) is 
an infinitesimal Hermitian matrix given by the flux of the covariant curl of 
A, 

V = Vx A- Ax A, (3.8) 

through the infinitesimal area bounded by C. (Note that because the com- 
ponents of A are matrices, the cross product A x A does not vanish, just 
as the cross product a x a = 2ia of Pauli matrices does not vanish.) This 
result does not extend to finite circuits, though, as there is no nonabelian 
version of Stokes' theorem. 

Like the ordinary geometric phase, U (C) may be interpreted as holon- 
omy, but now on a complex r-dimensional vector bundle rather than a line 
bundle. The base manifold is parameter space (M 3 in the examples consid- 
ered here), and the fibers are the degenerate eigenspaces spanned by |j(R)). 
The matrix-valued vector field A(R) is the connection form, which describes 
how states in the fibres are parallel transported along paths R(t) in the base. 
V(R) is the associated curvature form. 
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The nonabelian geometric phase does not depend on the choice of basis 
| j(R)) i n an y intrinsic way. Under a unitary change of basis 

r 

|j(R)>^XV jfc (R)|fc(R)>, (3.9) 
k=i 

the nonabelian version of the gauge transformation (2.8), U{C) transforms 
according to U{C) ->• WU(C)W^. This unitary conjugation represents sim- 
ply the change of basis applied to U(C) and leaves its observable properties 
unchanged. The fields A(R) and V(R) transform according to 

A -»• W AIU f - iVWW\ V -»• W VIU f (3.10) 

under (3.9). These transformation properties characterize the potential and 
field strength of Yang-Mills field theories, generalizing the electromagnetic 
analogy with the geometric phase (cf Section 2.2) to the nonabelian case. 

Given a particular closed circuit C, one can choose the basis |j(R)) 
so as to diagonalize U(C). Then under parallel-transport round C, each 
basis state returns to itself up to a phase factor. Under certain conditions, 
one can find a basis which diagonalizes every holonomy, regardless of C. 
In such cases, the degenerate and nondegenerate cases are essentially the 
same. In general, though, there is no choice of basis which renders every 
U{C) diagonal. What distinguishes the general case is that the generators 
of the symmetries responsible for the degeneracies depend on parameters 
themselves. A particular example, the A-doubling of electron doublets in 
diatomic molecules, is discussed by Moody et al. (1986). 

The discussion of geometric magnetism in the Born-Oppenheimer ap- 
proximation (cf Section 2.7) can also be generalized to the degenerate case. 
As before, we regard the parameters R as the coordinates of a slow system 
to which the original one, the fast system, is coupled. The Hamiltonian 
is given by P 2 /2M + H(R), where P is the slow momentum and the de- 
pendence on the fast coordinates and momenta is left implicit in H(R). 
We look for approximate eigenfunctions of the coupled system in which the 
fast system belongs to a degenerate multiplet of H(R). These are of the 
form \I/(r,R) = ^^ =1 ^(R)^' (r, R), where 0j(r,R) = (r|j(R)) are the 
fast eigenfunctions (which depend parametrically on the slow variables R). 
Substituting this ansatz into the stationary Schrodinger equation and mak- 
ing approximations similar to those in the nondegenerate case (the details 
are omitted), one sees that the slow system is described by a multicompo- 
nent wavefunction i/?(R) = (V'i(R), • • • , Yv(R)) which is an eigenfunction of 
the matrix-valued Hamiltonian 

H = (PI- A) 2 /2M + £(R)I. (3.11) 

Here A(R) is the matrix-valued vector potential of (3.5) and I is the identity 
matrix. Thus, coupling to a degenerate eigenstate introduces a generalized 
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form of geometric magnetism, in which the slow Hamiltonian acquires spin 
degrees of freedom which are coupled to a Yang-Mills-like potential. 

Finally, there is a connection between the nonabelian geometric phase 
higher-order degeneracies, ie points in parameter space where two or more 
degenerate levels coalesce. These points act as generalized monopoles sources 
for the nonabelian field V. 

3.3 Nonadiabatic evolution 

Geometric phases appear in contexts other than adiabatic evolution. An 
important generalization was found by Y Aharonov and J Anandan (1987), 
who showed that geometric phases accompany any cyclic evolution, not nec- 
essarily adiabatic, of a quantum state. Subsequently, J Samuel and R Bhan- 
dari (1988) showed that geometric phases can arise in evolutions which are 
nonunitary, not cyclic, and even discontinuous. What underlies these gen- 
eralizations is a natural prescription for the parallel transport of quantum 
phase factors. 

The setting for the Aharonov-Anandan phase is quite general. A state 
\ip(t)) evolves according to the Schrodinger equation (ihdt — H(t))\ij)) = 
under a (usually) time-dependent Hamiltonian H(t). After a time T the 
state is presumed to return to itself up to a phase factor, so that \ip(T)) = 
exp(i(3)\ip(0)). In the trivial case where H is time-independent and \ip) is an 
eigenstate, (3 is simply the dynamical phase —ET/h. Guided by this simple 
case, we are led to subtract from j3 a generalized dynamical phase, namely 
the time integral of the expectation value of energy 

Pd = ~\ J\m\^) dt = Im dt. (3.12) 

The phase which remains, f3 g d = n ft — fid, is given by 

P B = P-Pd = arg W0)|V(T)> - Im C dt (3.13) 

J o 

What makes the decomposition f3 = f3d + (3 g interesting is that f3 g is, in a 
sense to be explained, a geometric phase. First, (3 g depends only on the path 
through Hilbert space described by \ip(t)}, and not on the rate this path is 
traversed. Indeed, the state \4>(t)) = \ip(r(t))) describes the same path but, 
because of the reparameterization t — > r(t), at a different rate. However, it 
is easy to verify that its geometric phase (3 g , as given by (3.13), is the same as 
\ip(t)Ys, even though its dynamical phase is different (and it may no longer 
satisfy the Schrodinger equation). Second, /3 g does not depend on the overall 
phase of |Y>(i)). Indeed, if \ip(t)) is replaced by |Y>(i)) = exp(i0(t))|^(i)) in 
(3.13), a simple calculation shows that (3 g is unchanged. Such a substitution 
can be useful in evaluating (3.13). As \ip(t)) is periodic up to a phase, the 



23 



density matrix p^,(t) = \ip{t)){ij){t)\ is strictly periodic, and over the time 
< t < T describes a closed curve C in the space of density matrices. What 
the preceding discussion shows is that fi g is a geometric property of C. It 
depends not on the rate with which C is traversed, and it can be computed 
with (3.13), substituting any state \ip(t)) whose density \ijj(t))(ip(t)\ is equal 
to p^(t). 

The geometric phase for adiabatic cycles can be obtained as a spe- 
cial case of (3.13). Within the adiabatic approximation, a solution to the 
Schrodinger equation for a slowly cycled Hamiltonian H(R(t)) is given by 
| VC*)> = exp(i/3(t))|n(R(i))) ( cf ( 2 - 3 ))- From (3-12), the dynamical phase 
fid is given by —1/HJq E n (R(t))dt. The geometric phase is easily computed 
from (3.13) if \tp(t)) = |n(R(t))) is substituted for \tp(t)) (this is permissible, 
since the two differ by an overall phase). Since arg(n(R(0)|n(R(T))) = 0, 

fi g is given by — Im (ip\ip)dt = — Im <f>(n(R)\ Vn(R)) • dH, in agreement 
with (2.5). 

What underlies this general setting for the geometric phase is once more 
the mathematics of connections and fibre bundles. However, now all of 
Hilbert space is regarded as a fibre bundle! The base manifold % consists 
of all pure state density matrices p, normalized so that Trp = 1. (This 
is sometimes called the projective Hilbert space.) The fibre attached to p 
consists of all state vectors \ip) whose density matrices \ip){ip\ are equal to 
p. (Thus states in the same fibre differ only by a phase factor.) To an 
infinitesimal displacement dp there corresponds displacements \dip) in the 
fibre satisfying \ip+dip){ip+d r ip\ = p+dp. With a little algebra this reduces to 
\dip) — (ip\dip)\ip) = dp. This equation does not determine \dip) uniquely; to 
any solution can be added infinitesimal imaginary component ida\tp) along 
\ip). A connection fixes this component and determines a particular \dip) for 
a given dp. A natural choice is to require that 

(ip\dip) = 0. (3.14) 

This connection is precisely the condition (2.2) found by Berry for adiabatic 
cycles, but now generalized to arbitrary cycles in Hilbert space. Its holonomy 
is precisely the geometric phase fi g (3.13). 

This setting admits further generalizations. For example, one can calcu- 
late a geometric phase for a nonunitary evolution | </>(£)), for which {4>\4>) is no 

longer constant, simply by evaluating (3.13) with \ip(t)} d = n |0(t))/[(^|0)] 1 / 2 . 

In fact, one can calculate a geometric phase for an open path \ip{t)), 
where \ip(T)} and \ift(0)) are no longer proportional. The formula (3.13) for 
fig can be computed for open paths, and the result does not depend on the 
parameterization in time of \ip(t)), and remains unchanged by the change of 
phase \ip(t)) — > exp(i(j)(t))\tp(t)) . Further insight into the interpretation of 
the geometric phase for open paths is given in the following Section 3.4. 
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There is also a discrete version of the formula (3.13). Given a sequence of 
states |Y>i), |^2), • • • IV'w), one can associate a phase according to the formula 



Pg = arg(^i|^2)(V ; 2|V'3) • • • (V'aHV'i)- 



(3.15) 



One sees easily that (3 g does not depend on the overall phases of the states 
\tpj). If there are just two states, then f3 g = arg KY'ilV^)! = identically. 
But if N > 2, the phase j3 g is typically nonvanishing. Note that \iPn) need 
not be proportional to - the sequence of states need not close on itself. 

The formula (3.15) is the discrete analogue of (3.13), and defines a geo- 
metric phase for discontinuous evolutions. Thus, a geometric phase can be 
defined for processes which include measurements and wavefunction collapse. 
A version of formula (3.15) appears in work of Pancharatnam in optics, in 
which \ipj) represent polarization states. Panchartnam's contributions are 
discussed in Section 5.2. 

3.4 Geometric metric tensor 

How close are two quantum states \ip) and |</>)? The usual quantum mechan- 
ical measure of distance is the Hilbert space norm \ — (j)\ \ of their difference. 
There is another, first introduced by J Provost and G Vallee (1980), namely 
the quantity 1 — \ {il)\cj))\ 2 / {(f)\cj)) . Clearly this vanishes if the states are 
the same. What distinguishes it from the Hilbert space norm is that it is 
unchanged if \4>) or \4>) are multiplied by a complex scalar. In fact, what 
this quantity really describes is the distance between their normalized den- 
sity matrices p^ = / and p^ = \<p) {(j>\ / {(j)\4>) , and indeed can be 
alternatively expressed as 1 — Tr p^p^. For infinitesimally displaced states 
\<p) = \ip + dip), a little algebra gives 



This quantity defines a metric on the space of normalized density matrices 
H (or equivalently, the projective Hilbert space). 

The quantum metric enters into the analysis of the geometric phase in a 
number of interesting ways. Below we describe two aspects: a deeper under- 
standing of the geometric phase for open paths, and physically significant 
corrections to the Born-Oppenheimer approximation. 

The quantum metric determines geodesies, ie curves p(s) which mini- 
mize, or at least render stationary, the distance between two given density 
matrices in H . These are most easily described in terms of their correspond- 
ing state vectors; it may be shown that if \ip(s)) satisfies d 2 \ip)/ds 2 + \ip) = 
and if |^(0)), \ip'(0)) are normalized and orthogonal to one another, then the 
curve ptp(s) = \ip(s))(tp(s)\ is a geodesic with respect to the metric (3.16), 



ds 2 (dip) 



defn 

' (V#> <v#> 2 



(dip\I — p^\dtp) 



2 Tr ( d Pi>f ■ 
(3.16) 
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and the parameter s along the curve is its arclength. (The solutions \ip(s)} 
are of the form cos s \ip(0)) + sins ^'(O)).) 

A simple illustration is provided by two-state systems. In this case, H 
may be identified with the unit sphere (in the context of polarization optics, 
this is the Poincare sphere - see Section 5.2); each density matrix p may be 
expressed in the form (/ + e.er)/2, where e is a unit vector and cr are the 
Pauli matrices. The quantum metric (3.16) reduces to the usual metric on 
the sphere (up to a scalar factor), and the geodesies are the arcs of great 
circles. 

As shown in the previous Section 3.3, geometric phases can be defined 
via (3.13) for open as well as closed paths in %. The two cases may be re- 
lated as follows. An open path P may be closed by joining its endpoints with 
a geodesic. The resulting closed path, it turns out, has the same geometric 
phase as P. The geometric phase for continuous and discontinuous evolu- 
tions, (3.13) and (3.15), may be similarly related; given a closed sequence of 
density matrices pi, p2, ■ ■ ■ Pn , pi, we can construct a closed polygonal path 
p(s) in % by connecting its points with geodesies. p(s) turns out to have 
the same geometric phase as the discrete sequence pj. 

The quantum metric also plays a role in the Born-Oppenheimer descrip- 
tion of the dynamics of a slow system coupled to a fast one (cf Section 2.7). 
There it appears in the slow Hamiltonian as an additional scalar potential 
$ n (R) = ^2 a (ds 2 (\V a n(R)))/2M , which produces a 'geometric electric' 
force. The scalar potential appears at a higher order of adiabatic expansion 
than the vector potential A n (R); with appropriate scalings, the geometric 
magnetic force R x V n is of the order of the square root ^JrnjM of the mass 
ratio of the fast and slow systems, whereas the "geometric electric force" 
— is of the order m/M. However, the electric force can become signifi- 
cant in the vicinity of degeneracies R* (ie, points where the fast energy level 
E n (R) becomes degenerate). There the force is repulsive, and so tends to 
steer the slow dynamics away from points where the Born-Oppenheimer ap- 
proximation would break down. A discussion of the scalar potential can be 
found in the reviews of Berry (1989) and Jackiw (1988). Additional terms 
in the slow dynamics of comparable order to geometric electricity have been 
analyzed by Weigert and Littlejohn (1993). 

3.5 Beyond adiabaticity 

3.5.1 Corrections to the geometric phase 

For the adiabatically cycled eigenstates described by Eq. (2.3), the dynami- 
cal and geometric phases are but the first two terms in an asymptotic expan- 
sion for the total phase. An elegant iterative scheme for evaluating higher 
order corrections was developed by MV Berry (1987c). These corrections to 
the geometric phase turn out to be geometric phases themselves. The full 
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asymptotic series diverges in a manner which reveals in a fundamental way 
the breakdown of the adiabatic approximation. 

The method is illustrated by the example of a spin- 1/2 particle in a 
magnetic field (cf Section 2.3). We consider a time-dependent Hamilto- 
nian Hq(t) = ^ / uBo(r) • <x, cyclic over the infinite interval [—00,00] with 
Bo(±oo) = B*z. Bo(t) is analytic, and we assume that Bq(t) > B m i n 
for all r, so that the Hamiltonian does not pass through any degeneracies, 
r = et is a dimensionless time, and e a small parameter (with dimensions 
of frequency) describing the slow rate of change of the Hamiltonian. (The 
corresponding dimensionless small parameter is 5 = eh/ fiB m i Q .) 

Let us consider the evolution of a spin state \ip(t)} with initial condition 
00)) = |— ). The analysis is facilitated by transforming to a rotating 
frame in which the direction £>o of the magnetic field is fixed (say, for def- 
initeness, along z). Such a frame is described by a rotation matrix 1Zq{t) 
satisfying 

bo(r) = fto(r) ■ z. (3.17) 

This condition does not determine 1Zq{t) uniquely, but only up to a rotation 
(applied from the left) about bo- Differentiating with respect to r (this is 
denoted by a dot ') , we obtain a condition 

6 (r) = u> (r) x b (r) (3.18) 

on the angular velocity cjo of the frame, where ojqi = \ X^jfci e ijkT^okiT^oji 
(eijk is the Levi-Cevita symbol). (3.18) does not determine uq uniquely, 
but only up to its component along bo- We fix this indeterminacy (and the 
corresponding indeterminacy in TZq) by requiring that 

u> -b = 0. (3.19) 

(3.19) is a condition for parallel transport, analogous to the requirement of 
moving your arm without unnecessarily twisting your wrist (cf Section 3.1). 

In the rotating frame, the Hamiltonian H\{t) is given by ^/zBi(t) • er, 
where 

B 1 (t) = B (t)z-B l1 {t) (3.20) 

is the effective magnetic field. B^i(r) = £/i/2^7^q" 1 (t) -uq(t) is the Larmor 
field induced in the rotating frame (the magnetic analogue of the Coriolis 
force). As the rotation is slow, the Larmor field is small (of order e), and in 
the lowest order adiabatic approximation (cf Sections 2.1, 2.2) it is simply 
neglected. Then H\(t) is diagonal, and the time-dependent Schrodinger 
equation is easily solved. The geometric phase 7_(Co), half the solid angle 
of the cycle Co described by Bo(r), appears as a direct consequence of the 
parallel transport condition (3.19) when the solution is transformed back to 
the fixed frame. 
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If the Larmor field is not neglected, then Bi(r) may be regarded as a 
slowly- varying magnetic field to which the above procedure may be applied 
as before. That is, we transform to a second rotating frame TZ\{t) in which 
bi(r) is fixed along z. The frame is uniquely determined by imposing the 
parallel transport condition u)\ ■ bi = 0. Note that because hi differs from 
z by terms of order e, the angular velocity uj\ is of order e. In this second 
frame, H2(t) = ^B2(r) • <r, where the effective magnetic field B2(t) is 
given by B\o z — B^ 2 (t), with B^ 2 (t) = eh/2[iTZ~ 1 (T) ■ u>\(t). The new 
Larmor field B^2 is of order e 2 . 

It is easy to see that this procedure can be performed iteratively; we 
can transform through successive rotating frames TZ u (t) to compensate the 
Larmor fields Bi„(r) induced in the preceding ones. The geometric phase 
and its corrections emerge when the sequence of rotations is unravelled and 
the solution transformed back to the fixed frame. One obtains a formal 
series 

oo 

7-(C ) + X>-(C fc ), (3-21) 
k=i 

in which the corrections to the geometric phase are themselves geometric 
phases; 7_(Cfc) is half the solid angle of the cycle described by B^(r). 

It turns out that the series (3.21) is divergent. Although successive 
cycles (7^ initially decrease in extent by a factor of e, asymptotically the 
terms 7_(Cfc) behave as e k k\, a form characteristic of divergent asymptotic 
series (Dingle 1973). An optimal approximation is obtained by truncating 
the series at its smallest term k ~ 1/e, for which e k k\ ~ exp(— 1/e) is 
exponentially small. 

That the series (3.21) diverges can be understood on general grounds. 
If, on the contrary, the series converged, along with the infinite sequence of 
rotating frames, than upon transforming back to the fixed frame, we would 
obtain an exact solution \tp(t)) of the Schrodinger equation which would be- 
gin and end, at t = — oo and t = +oo respectively, in the spin-down state 
| — ). This would imply a vanishing transition probability (+|VK+oo)). But it 
is known that, with the exception of some special Hamiltonians, the transi- 
tion probability is not zero, but instead (for Bo (t) analytic) is exponentially 
small, of order exp(— 1/e) (Dykhne 1962). These exponentially small transi- 
tions (the subject of the following section) signal the ultimate failure of the 
adiabatic approximation, and cannot be captured by the iterative scheme. 
It follows that this scheme must inevitably diverge. 

3.5.2 The geometric amplitude 

Adiabatic evolution is an approximation. For typical slowly varying Hamil- 
tonians H(t), where r = et, no matter how small the rate of change e, 
there are transitions between evolving eigenstates. The transition probabil- 
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ity vanishes as e goes to zero, and if H depends analytically on r it does so 
exponentially rapidly (faster than any power of e). 

In a study of adiabatic evolution in two-state systems, MV Berry (1990b) 
found that this transition probability may contain a geometric factor. The 
setting is similar to that of the preceding Section 3.5.1. We consider a 
two-state Hamiltonian H(r) = R(t) • <r, where R(r) is analytic in the di- 
mensionless time r = et and is nonvanishing for real r. As r — > ±00, R(t) 
approaches limiting values R±* which need not coincide, so that R(r) need 
not be closed. From (2.3), the adiabatic solution for the up-state is given 
by \ip(t)) = exp(i<5+(i) + i^ + (t))\ + (R(r))), where the eigenstates | ± (R)) 
are given by (2.14) and the evolving dynamical and geometric phases 5 + (t) 
and 7+(i) are given explicitly by 



S+(t) = — - R(r)dT, 7+(t) = — / (1- cos 6) <j) dr. (3.22) 



Here (R, 9, 4>) denote the polar coordinates of R(t), and the dot ' denotes the 
derivative with respect to r. We have used (2.15) for the vector potential 



Nonadiabatic transitions may be attributed to roots of R(t) in the com- 
plex r-plane. The dominant contribution to the transition amplitude comes 
from the root r c closest to the real axis. (We assume r c is a simple zero of 
R 2 {t), so that it is a square-root branch point of R{t).) Davis and Pechukas 
(1977) showed that the transition amplitude may be obtained by analyti- 
cally continuing the adiabatic solution around the branch point r c . (Under 
this continuation, the up- and down-states | + (R(r))) and | — (R(r))) are in- 
terchanged.) The transition probability = |(— (R(oo))|'0(oo))| 2 is given 

by a product of factors exp(— | A + |) exp(r + ), whose exponents are obtained 
from the analytical continuation of the dynamical and geometrical phases. 
Explicitly, we have that 



In the expression for r+, the contour is taken round the cut joining r c and 
its complex conjugate r*. 

The factor exp(r + /2) is the geometric amplitude. It depends only on 
the cycle followed by R(t), not on the slowness parameter e nor on h. It 
can be shown to vanish if R(r) lies in a plane, or if R(t) can be rotated into 
itself about an axis through the origin. Thus it vanishes for the well-known 
Landau-Zener model R(r) = (a, a, At) (Zener 1930), for which R(r) de- 
scribes a vertical line. A simple example with nonzero geometric amplitude 
is R(r) = (acos(ojr 2 ), a sin(wr 2 ), At). It describes a helix on the cylinder 
of radius a which winds for r < and unwinds for r > 0. For this curve r + 
is given by — a 2 wsgn A /A 2 . 





A+(R). 




(3.23) 
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The geometric amplitude was observed by Zwanziger et al (1991) in a 
nuclear magnetic resonance experiment (cf Section 6.2). A time-dependent 
spin Hamiltonian with the requisite properties was generated by applying a 
modulated rf magnetic field to a sample of spin-1/2 nuclei, and the induced 
magnetization was observed in a rotating frame. The logarithm of the tran- 
sition probability was plotted as a function of the slowness parameter 

e. Theory predicts that this plot should be linear (due to the dynamical 
amplitude) with a nonzero intercept (due to the geometric amplitude) . This 
behaviour was confirmed in the experiment, with good quantitative agree- 
ment between the observed and predicted value of T + . 

4 Classical Systems 

Following Berry's discovery of the geometric phase, JH Hannay (1985) found 
an analogous phenomenon in classical mechanics. Adiabatically cycled clas- 
sical oscillators undergo a shift in their phase of oscillation. Part of this 
shift is simply the time integral of the frequency. But there appears in ad- 
dition a purely geometrical contribution, now known as the Hannay angle. 
The Hannay angle includes and generalizes such well-known phenomena as 
the precession of the Foucault pendulum. In many cases it may be directly 
observed. The Hannay angle emerges as the classical (h — > 0) limit of the 
Berry phase (in cases where chaos is absent from the classical dynamics). 
In the context of coupled systems, it leads to classical geometric magnetism 
in the slow dynamics. 

F Wilczek and A Shapere (1989b-1989d) studied phenomena in classical 
kinematics similar to the Hannay angle. A simple example is the manner 
in which falling cats and platform divers change their spatial orientations 
through a cycle of internal motions. Provided their (conserved) angular mo- 
mentum is zero, the change in their orientation does not depend on the rate 
at which the internal motions are performed. Classical geometric phases 
of this type have been investigated systematically by J Marsden, R Mont- 
gomery and coworkers (Marsden et al. 1991, Marsden 1992) as part of a 
general theory of Hamiltonian dynamics with symmetry. 

4.1 The Hannay Angle 

4.1.1 The Classical Adiabatic Theorem 

Imagine a swinging pendulum slowly shortened to half its original length. 
During this process, is anything conserved, either exactly or approximately? 
This question is addressed by the classical adiabatic theorem. While energy 
is not conserved (work is being done), there is an approximately conserved 
quantity, the classical adiabatic invariant, the action. 
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The classical adiabatic theorem is most simply described for systems 
of one degree of freedom, so we concentrate on this case. The dynamics 
in the two-dimensional phase plane is described by Hamilton's equations 
q = d p H,p = —d q H. The classical Hamiltonian H(q,p,H) is taken to 
depend on parameters R (like the quantum Hamiltonians of Section 2). If 
these are held fixed, then energy is conserved, and orbits lie on contours 
of H. Let us assume the motion is bounded, so that the energy contours 
are closed. The orbits go round them periodically with oscillation frequency 
w = {9eI)~ 1 -, where 



is the action. Apart from the factor of 2ir, it is just the area enclosed by the 
orbit of energy E. If the parameters vary in time, then energy is no longer 
conserved. However, if the variation is slow compared to the oscillation 
frequency uj, then the action (4.1) is nearly conserved, and the variation of 
the energy is determined implicitly by 



This is the classical adiabatic theorem. 

In higher dimensions, the classical adiabatic theorem generalizes pro- 
vided the dynamics is either integrable or ergodic. Integrable systems are 
characterized by a high degree of symmetry. They possess N independent 
constants of the motion (where N is the number of degrees of freedom) which 
are in involution (that is, the Poisson bracket between any pair of them van- 
ishes). Typical orbits in an integrable systems execute periodic or quasiperi- 
odic motion on iV-dimensional tori embedded in the 2iV-dimensional phase 
space, and the equations of motion may be solved more or less explicitly. 
In sharp contrast stand the ergodic systems, for which energy is the only 
constant of the motion. Typical orbits cover the (2N — l)-dimensional en- 
ergy surface uniformly, so that time averages and microcanonical ensemble 
averages coincide. The adiabatic theorems in higher dimensions are rather 
more difficult to state and considerably harder to prove. The difficulty is 
related to the fact that the adiabatic invariants are not as well conserved as 
in the one-dimensional case (Lochak and Meunier 1988). 

4.1.2 Basic account of the Hannay angle 

The classical adiabatic theorem (4.2) determines the time evolution of the 
energy of a one-dimensional oscillator under a slowly changing Hamiltonian. 
J Hannay (1985), motivated by Berry's analysis of the quantum geometric 
phase, asked a more refined question: what is the evolution of the phase of 
oscillator? 




dq dp 



(4.1) 



I{E{t)Mt))=h = /(£(0),R(0)). 



(4.2) 
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To frame this question precisely we introduce canonical coordinates, the 
action-angle variables {9,1). These are related to q,p by a canonical trans- 
formation q = q(9,I, R),p = p(0, 1, R) which depends on parameters. The 
action /, the generalized momentum, determines the energy of the orbit 
according to (4.1). The angle 6, the generalized coordinate, describes the 
phase of the oscillator, and varies periodically around the orbit with period 
2-7T. Expressed in terms of action-angle variables, the Hamiltonian H(I, R) is 
independent of 9, so the canonical equations of motion, 9 = djH, I = —dgH, 

imply that / is constant while 6 evolves at a constant frequency oj = n djH. 

Now suppose the parameters R of the Hamiltonian are taken slowly 
round a closed cycle C, returning to their original values R(0) after a time 
T. Since the action is nearly conserved, a trajectory returns approximately 
to the contour of H (I, R(0)) that it started on. Of course, it will be displaced 
by an angle A9 along the orbit from its initial position. Hannay obtained 
the two dominant contributions to AO (the corrections are of order 1/T). 
The first is the dynamical angle 6d = f uj(I ,~R(t))dt. It is simply the 
time integral of the frequency, taking into account its slow variation with 
parameters. The second, unexpected, contribution, 



is the geometric angle, or the Hannay angle. It is expressed in terms of the 
line integral of a vector field A(I, R) given by 

A(J, R) = i- j** p(0, I, R) Vg(0, 1, R) dO. (4.4) 

The Hannay angle is geometric in the same sense as the geometric phase; it 
depends only on the parameter cycle C, and not on the rate it is traversed. 
Using Stokes' theorem, the Hannay angle can also be expressed as the surface 
integral —dj f s Y(I, R) • dS, where S is a surface bounded by C and V = 
V x A. From (4.4), 



As with the geometric phase, the presence of time-reversal symmetry (cf 
Section 2.5) imposes constraints on the Hannay angle. In particular, if the 
orbits themselves are time-reversal symmetric (for all values of parameters) , 
then the Hannay angle vanishes identically. 

Like the geometric phase, the Hannay angle may be interpreted as the 
holonomy of a connection on a fibre bundle (cf Section 3.1.3). The param- 
eters R constitute the base manifold, and the fibres attached to them are 
orbits of the Hamiltonian H(q,p, R) with given action /. These orbits are 




(4.3) 




(4.5) 
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parameterized by the angle 9, and the connection <9/A(7, R) describes how 
9 changes under parallel transport. 

The Hannay angle can be detected by determining the initial and final 
angles of an adiabatically cycled oscillator, and subtracting from their dif- 
ference the dynamical angle 9^. As 9d depends only on the parameter cycle 
R(i), it can be computed without solving the equations of motion explic- 
itly. In this way, Hannay's analysis, together with the adiabatic theorem 
(4.2), provides a complete asymptotic solution of the equations of motion 
for slowly varying Hamiltonians in one dimension. The analysis generalizes 
directly to integrable systems in higher dimensions, where there are as many 
Hannay angles as there are degrees of freedom. However, detecting them can 
be a more delicate matter, as the adiabatic theorem and corrections thereto 
become more complicated in higher dimensions (Golin et al. 1989). For er- 
godic systems, there is no general analogue of the Hannay angle, though a 
related holonomy phenomenon emerges in some special cases (Robbins 1994, 
Jarzynski 1995, Schroer 1996). 



4.1.3 Examples 

A bead moving freely round a planar wire loop provides a simple example 
of the Hannay angle. Suppose the bead has unit mass, and the loop has 
circumference C, area A, and is otherwise arbitrary in shape. While the 
system as a whole has time-reversal symmetry, the individual orbits do not 
(they are either clockwise or counterclockwise). The action I is just the 
velocity v, and the angle 9 = 2ns /C is proportional to the distance s along 
the wire. Suppose the loop is slowly rotated once around an interior point 
over a time T, and let As denote the total distance the bead has travelled 
during this time. The dynamical contribution Sd (which is C/2ir times the 
dynamical angle 9d) is just vT (the adiabatic theorem implies v is nearly 
constant). The geometric contribution Sd may be computed to be — AttA/C 2 . 
It accounts for the distance the bead has slipped as the wire turned. For a 
fixed circumference C, Sd assumes its maximal value, C, for a circular loop; 
in this case, the bead simply slips a full revolution. For thin, elongated 
loops, Sd approaches zero; the only slippage occurs at the hairpin turns. 

The second example is the family of harmonic oscillators H = [Z — 
X)p 2 /2 + Ypq + (Z + X)q 2 , where the term in qp breaks time-reversal sym- 
metry. For X 2 + Y 2 < Z 2 , the motion is bounded and oscillatory, and we 
consider adiabatic cycles confined to this region of parameter space. The 
flux field V of (4.5) may be computed to be 

V(/ ' R) = 1 (^-Xyy (46) 

Note that because V is linear in /, the Hannay angle 9 g is independent of 
/. 9 g may be described geometrically as follows. We regard the parameter 
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space as three-dimensional Minkowski space, in which X and Y are space- 
like coordinates, Z is time-like, and the region X 2 + Y 2 < Z 2 is the interior 
of the light-cone. Then g (C) is the area of the projection of the closed 
curve C onto the unit mass shell Z 2 — X 2 — Y 2 = 1. 

The precession of a Foucault pendulum is yet another example of a Han- 
nay angle. Consider a spherical pendulum initially set to oscillate in a plane. 
During the course of a day, the axis of gravity (whose direction represents 
the parameters of the pendulum) is carried once round a line of latitude 
a. At the day's end, the plane of oscillation has turned through an angle 
27r(l— cos a). This is just the solid angle described by the axis of gravity, and 
may be computed as a Hannay angle using (4.3). In fact, the phenomenon 
described by the Hannay angle is more general. If one could adiabatically 
but otherwise arbitrarily turn the earth around, returning it to its original 
orientation, the plane of oscillation of the spherical pendulum would still 
precess through an angle equal to the solid angle subtended by the earth's 
axis. The spherical pendulum is a two-degree of freedom system, integrable 
because of its axial symmetry. Therefore, it has two Hannay angles, one 
each for its polar and azimuthal degrees of freedom. For planar oscillations, 
however, only the azimuthal Hannay angle, which describes the Foucault 
precession, is nonzero. 

Let us briefly mention two final examples drawn from plasma physics 
and astronomy, respectively. Littlejohn (1988) found a Hannay angle in the 
cyclotron oscillations of charged particles in strong magnetic fields. Berry 
and Morgan (1996) (see also Golin et al. (1989)) compute "the Hannay 
angle of the world" , an annual displacement of the earth along its orbit of 
approximately —150 meters (or —2 x 10~ 4 seconds of arc) caused by the 
adiabatic and periodic variation of Jupiter's gravitational field. 

4.1.4 Semiclassical limit 

Following Hannay's work, Berry (1985) showed that the Hannay angle emerges 
in the semiclassical limit of the geometric phase. The precise relation be- 
tween the two is given by 

7n + i(c) - 7b (o ^ ; n 7°°) e g (i,c), (4.7) 

nn=l 

in which classical action is related to the quantum number by the semiclas- 
sical quantization condition / = nh. (4.7) conforms to the usual relation 
between quantum phases and classical angles, a more familiar instance be- 
ing the relation (E n+ \ — E n )/h — > uj between energy levels and frequencies, 
whose time integral connects the dynamical phase and the dynamical angle. 
The result (4.7) is derived by substituting the WBK approximation for the 
eigenstates |n(R)) into the formula (2.4) for the quantum vector potential 
A n . One obtains thereby the expression (4.4) for the classical vector poten- 
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tial. This analysis extends to integrable systems in more than one dimension, 
for which semiclassical eigenfunctions may be similarly determined. 

On the other hand, for systems whose classical dynamics is chaotic, the 
quantum geometric phase has no classical limit in general. However, the 
vector field V n , when averaged over quantum number, does have a classical 
limit (Robbins and Berry 1992), which turns out to play a role in the classical 
Born-Oppenheimer approximation (cf Section 4.1.5). The classical limit 

V(E,R) = - 2 g E ^ E R) d E (d E n(E,-R)J™{VH t (q, P ,-R) x VH(q,p,K))dt 

(4-8) 

is expressed in terms of the time integral of an antisymmetric correla- 
tion function of the parameter gradient of the Hamiltonian, Vff(g,p,R). 
Q(E, R) = fn( qp n-)>Ed N qd N p denotes the phase volume inside the energy 
shell H = E, {■) denotes the normalized microcanonical average on the en- 
ergy shell, and Vi7< denotes the function Vi7 evolved for a time t with the 
classical dynamics. In one dimension, (4.8) is equivalent to (4.5). 



4.1.5 Classical reaction forces 

As in Section 2.7, we suppose the parameters R are no longer externally 
prescribed, but instead are the coordinates of a slow system to which the 
original one, the one-dimensional oscillator of Section 4.1.1, is coupled. Both 
systems are governed by classical mechanics. We take their Hamiltonian to 
be P 2 /2M + H(q,p, R), where P is the slow momentum. The Hannay angle 
describes a holonomy in the fast dynamics induced by the slow dynamics. 
As in the quantum case, there is a reciprocal effect on the slow dynamics 
induced by the fast dynamics. This effect is found by averaging over the fast 
motion. According to the adiabatic theorem, the action of the fast system 
is nearly constant. In the lowest order approximation, we treat the action 
as strictly constant. Then the energy H(I,~R) acts as an potential in the 
slow dynamics, whose effective Hamiltonian is given by P 2 /2M + H(I, R). 
This is the classical analogue of the Born-Oppenheimer approximation. 

First-order corrections account for the fluctuations of the action about 
its mean. These generate an effective vector potential in the slow Hamilto- 
nian, so that P — > P — A(I, R). The vector potential is precisely the field 
(4.4) whose line integral gives the Hannay angle. This is classical geometric 
magnetism. 

One example of geometric magnetism is the precession of a symmetric top 
(Berry and Robbins 1993a). In this case, the top's axis is the slow system, 
and the rotation about this axis is the fast system. Another example from 
plasma physics is the drift of the guiding centre of electrons in strong, slowly- 
varying magnetic fields (Littlejohn 1988). The guiding centre, ie the centre 
of the instantaneous Larmor orbit, constitutes the slow degrees of freedom 
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(it remains fixed in constant fields), and the motion about the guiding centre, 
the fast degree of freedom. 

Geometric magnetism generalizes to higher dimensional fast systems, 
both integrable and chaotic. In the latter case, geometric magnetism is 
partnered by an additional dissipative force, deterministic friction; the slow 
dynamics is no longer conservative, as the chaotic motion of the fast system 
acts as an effective heat bath. The magnetic and dissipative forces appear 
as the antisymmetric and symmetric parts of a linear response tensor (Berry 
and Robbins 1993b, Jarzynski 1993). 

4.2 Symmetry and kinematic geometric phases 

Kinematic geometric phases arise in systems with symmetry. They describe 
the evolution in degrees of freedom which the symmetry acts upon (eg, 
orientational degrees of freedom in case of rotational symmetry). Their 
structure depends only on how the symmetry is realized, and not on the 
details of the dynamics. The geometric character of their evolution is exactly 
enforced by the equations of motion, and does not depend on assumptions of 
adiabaticity. In the following sections, two different treatments of kinematic 
geometric phases are described. 

4.2.1 Falling cats and swimming paramecia 

Even with nothing to push against, bodies in free fall, such as cats and plat- 
form divers, can change their orientation through internal motions. Perhaps 
counterintuitive at first, this phenomenon is a direct conseqence of conser- 
vation of angular momentum, and can be calculated as such. The main 
difficulty in doing so is in the definition of the body's orientation. This 
problem does not arise for truly rigid bodies, whose every configuration is 
related to a chosen reference by a unique rotation. But for an body which 
can bend and twist and stretch, ie which has internal degrees of freedom, it 
is not clear which configurations are "upright" and which ones are "tilted". 

There are various prescriptions for determining orientation. One is to re- 
fer to the orientation of the principal axes of inertia. Another is through the 
Eckart frame, a common convention in molecular physics. Both presciptions 
have particular advantages; neither is well defined in all cases. 

In fact, the structure of the solution emerges most clearly if the definition 
of orientation is made more or less arbitrarily. This is the point of view taken 
by Shapere and Wilczek(1989b). For simplicity, suppose the body consists 
of a finite number N of masses m a at positions r a . No assumption is made 
concerning their dynamics, except that they obey Newton's laws and that 
there are no external forces. By choosing a suitable reference frame, we may 
assume that the centre of mass ^ a m a r a /^ a m a begins and remains at 
the origin. 
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By some means, certain configurations (Ri, . . . , Rat) are designated to 
be "upright". This designation must be consistent, in that if two configu- 
rations are related by a rotation, then only one of them can be upright. It 
must also be complete, so that every configuration can be made upright by 
applying a rotation. If both these conditions hold, then every configuration 
(ri, . . . ,r n ) can be obtained by applying a unique rotation 1Z to a unique 
upright configuration (Ri, . . . , Rat). The body's motion, 

r a (t)=K(t)--R a (t), a = l,...,N (4.9) 

is then determined by the evolution of its orientation TZ(t) and its internal 
dynamics Ri(i), . . . , Rat(£)- In the language of rigid body dynamics, R a (i) 
describes the motion in the body frame, and r a (t) the motion in the space 
frame. 

An object fixed in the space frame appears to rotate with angular velocity 
— fl(t) in the body frame. fl(t) is defined by the relation 

n-R = n- (ft x r), (4.10) 

which can be made to hold for arbitrary vectors R. 

Regardless of the details of the dynamics, the total angular momentum 
1 = J2 a m oJa x v a in the space frame must be conserved. Conservation of 
angular momentum leads to the relation 

nit) = i' K (^-Yl m « R « x ( 4 - n ) 

between the angular velocity fi, angular momentum L = ■ 1, velocities 
V„ = R Q , and inertia tensor Iy = J2a m a(Ra^ij ~ RaiRaj) in the body 
frame. Given the internal dynamics R a (t), the relations (4.10) and (4.11) 
may be used to determine the orientation lZ{t). 

The case of zero angular momentum (1 = L = 0) is particularly inter- 
esting. Then one can show that the orientation TZ(t) does not depend on 
the rate of the internal dynamics. (One can, in effect, multiply through 
by dt in (4.10) and (4.11).) In particular, if the internal motion is peri- 
odic along a cycle C, then after one period it generates an overall rotation 
TZg = 7£(T)7£(0) -1 of the body. The rotation 1Z g depends only on the cycle 
C, and not on the rate it is traversed. It is a purely geometrical quantity, a 
nonabelian analogue (rotations do not commute) of the Hannay angle. It is 
the means by which a cat dropped upside-down manages to land on its feet. 

Like the other geometric phases considered so far, this can be described 
in terms of holonomy on a fibre bundle. The fibre bundle is the configuration 
space of the masses m a , constrained to have centre of mass at the origin. 
The base manifold consists of the distinct internal configurations, no two 
of which can be made to coincide by a rotation. (Mathematically, this is 
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the quotient of the configuration space by the rotation group.) It can be 
parameterized by the upright configurations, provided these are always well- 
defined. The fibre consists of the different possible orientations of a given 
internal configuration, and may be parameterized by the rotation group. 
A connection determines an infinitesimal rotation 51Z for an infinitesimal 
change of internal configuration 5R a , and the physically relevant connection 
is determined by conservation of angular momentum. Parallel transport 
round a closed cycle C in the base yields the holonomy lZ g . The curvature 
of this connection, which describes the net rotation after an infinitesimal 
cyclic change of shape, may be computed using a formula analogous to 
(3.8). 

The designation of upright configurations is arbitrary, and one can always 
choose them differently. It can be verified that the holonomies 1Z 9 do not 
depend intrinsically on this choice, and that the various quantities which 
describe the orientation transform sensibly when a different choice is made. 
In particular, the curvature and connection transform according to equations 
similar to (3.10). 

A more complex example studied by Shapere and Wilczek(1989c) con- 
cerns swimming through viscous fluids at low Reynolds number. High vis- 
cosity and low velocity means that bodies remain at rest unless acted upon 
by a force - a realization of Aristotelian mechanics! Bodies locomote by 
executing a series of internal motions, ie swimming strokes. Conservation of 
linear and angular momentum determine the consequent changes in position 
and orientation, but now the motion of the fluid must be taken into account. 
Solving the fluid dynamical problem leads to a relation similar to (4.11) for 
both the linear and angular velocity in terms of rates of change of shape. 
Applications in biology include locomotion of micr organisms. 

4.2.2 Geometric phases and reduction 

Kinematic geometric phases can be framed in terms of the general theory of 
symmetry in Hamiltonian systems, called reduction. This formulation, ex- 
tensively developed by J Marsden, R Montgomery and coworkers (Marsden 
et al. 1991, Marsden 1992), has led to a unified perspective on a variety of 
geometric phase phenomena as well as a number of interesting examples. 

Dynamics is simplified by the presence of symmetry, and Hamiltonian 
dyanamics especially so. A symmetry always reduces the dimension of a 
dynamical system by one, but in the Hamiltonian case, it reduces it by two. 
Indeed, when there is sufficient symmetry, the systems are integrable, and 
the dynamics can be solved analytically (cf Section 4.1.1). 

The simplest symmetry is translation invariance in an ignorable coor- 
dinate qn ■ Hamilton's equation p^ = —dH/dqw = (more generally, 
Noether's theorem) implies that the conjugate momentum pat is conserved. 
Fixing its value pat = pno leads to a Hamiltonian H(qi, . . . , qjy-i,Pi, ■ ■ ■ ,pn-i, 
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in two fewer variables. Once their dynamics is obtained, the evolution of the 
ignorable coordinate may be determined by integrating Hamilton's equation 
Qn = d PN H(q 1 , q N -i,Pi, ■ ■ ■ ,Pu-i,Pno)- The theory of reduction ex- 
tends this analysis to general symmetries. As in the simple example, the 
analysis consists of two stages. The first is reduction, wherein the conserved 
momenta are fixed, the ignorable coordinates ignored, and Hamilton's equa- 
tions expressed in terms of a reduced set of coordinates. In all but the 
simplest cases, finding these coordinates is a nontrivial task (Jacobi's "elim- 
ination of the node" in the rotationally symmetric iV-body problem is a 
representative example.) The second state is reconstruction, wherein the 
evolution of the ignorable coordinates is determined from the reduced dy- 
namics. It is in the reconstruction stage than kinematic geometric phases can 
appear. It turns out that the equations of motion for the ignorable coordi- 
nates can be separated into dynamical and geometric terms. The geometric 
term depends only on the path described by the reduced dynamics, and not 
on the rate it is traversed. The dynamical term, by constrast, does depend 
on the rate. Moreover, the form of the geometrical term depends only the 
realization of the symmetry, and not on the details of the Hamiltonian. 

In the example discussed in the preceding section, the flexible body com- 
posed of N point masses, the ignorable coordinates correspond to 1Z, the 
body's orientation, the conserved momentum to the space-fixed angular mo- 
mentum 1, and the internal coordinates R Q (i) to the reduced system. The 
evolution of R(i) is obtained by integrating the linear system (4.10), in 
which angular velocity ft(t) is given in terms of the reduced dynamics by 
(4.11). Eq. (4.11) is an example of a reconstruction formula divided into 
its dynamical and geometric parts. The first term, J -1 • L, is the dynamic 
contribution. The second, I -1 • (%2 a m Q R a x V a ), is the geometric contribu- 
tion. As discussed above, its integral depends only on the path C described 
by the internal dynamics R Q (t). 

R Montgomery (1991) found an even simpler example in the dynamics 
of a free rigid body. Its configuration (with its centre of mass fixed at 
the origin) is described by a rotation 1Z, for which it will be convenient to 
use the Euler angle parameterization, 1Z = TZ z (4>)TZ y (9)Tl z (if;). The space- 
fixed angular momentum 1 is a constant of the motion, and without loss of 
generality we can assume it lies in the z direction. The body-fixed angular 
momentum L(t) = TZ(t) • 1 has polar coordinates (6(t),4>(t)) on the sphere 
of constant radius radius L = |L| = |1|. Its dynamics is determined by the 
Euler equations L = x L, where fi = I^ 1 • L is the angular velocity and / 
is the inertia tensor. The Euler equations describe the reduced dynamics, in 
which L, or equivalently 6 and 4>, are the reduced coordinates. The ignorable 
coordinate is ip. 

Under the Euler equation, L executes a periodic motion with period T. 
The fact that L(T) = L(0) implies that rigid body returns after one period 
to its original orientation, up to a rotation about the direction L(T). Indeed, 
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one has that K(T) = K(0)K z (Aip) = ^ L(T) (A^)7e(0), so that the amount 
of rotation, Aip = ip{T) — il>(0), is just the shift in the ignorable coordinate 
ip. With a little algebra, the equation of motion for ip may be written in the 
form 

ip = 2E + cosO<j), (4.12) 

where E = |L • I^ 1 ■ L is the conserved rotational energy. (4.12) is another 
example of a reconstruction formula divided into dynamical and geomet- 
ric contributions. In this case, the dynamical contribution integrates to 
2ET, where as the geometric contribution gives (modulo 2ir) the solid angle 
§ c cos 6 dip described by the body-fixed angular momentum L(i) over a single 
period. 

Another example, due to Alber and Marsden (1992), concerns phase 
shifts in soliton collisions. A spectacular property of certain nonlinear partial 
differential equations, such as the Korteweg-de Vries equation Ut + 6UU X + 
U X xx = 0, is the existence of solitons. These are waveforms U(x — vt) which 
propagate at a fixed velocity v whilst maintaining their shape, the nonlinear- 
ity in the wave dynamics just balancing the dispersion. When two solitons 
collide, they emerge after a complicated but transient period of interaction 
essentially intact. The only signature of their collision is a shift Ax in their 
position x ± Ax — vt (+ for the faster soliton, — for the slower one). This 
shift, which may be determined from an asymptotic analysis of the exact 
solutions, can be interpreted as a geometric phase. The same analysis ap- 
plies to a variety of other integrable nonlinear partial differetial equations, 
including the nonlinear Schrodinger equation and the Klein-Gordan equa- 
tion. 

4.2.3 Optimization and control 

The study of kinematic geometric phases suggests a problem in optimization. 
What is the most efficient internal motion for producing a given holonomy? 
As a specific example, consider a satellite orbiting the earth in a fixed ori- 
entation, so that its angular momentum about its centre of mass is zero. It 
is desired to change its orientation using as little energy as possible. The 
obvious method of applying external torques (supplied by jets attached to 
satellite, for example) has the drawback that any uncompensated force or 
torque leaves the satellite in a state of translation or rotation. The alterna- 
tive is execute a sequence of internal motions (eg, flywheels driven by motors, 
for example); in the manner described above, these generate a change the 
orientation without altering the total linear or angular momentum. 

The problem of finding the most efficient cycle with a given holonomy 
leads to a constrained variational problem in the internal degrees of freedom. 
It may also be formulated in terms of control theory. The solution depends 
very much on how efficiency is measured. In certain circumstances, it may 
be advantageous to make some sizeable internal motions initially in order 
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to achieve a configuration for which the geometric phase curvature form (cf 
Section 3.1.2) is large; then much smaller subsequent motions can produce 
relatively large effects. In other situations, it may be necessary to restrict 
to small motions about some equilibrium configuration; in this case, the 
variational problem may be analyzed in terms of infinitesimal circuits. The 
solutions often have a physical interpretation in terms of motion in a mag- 
netic, or more generally, a Yang-Mills field. Shapere and Wilczek (1989d) 
find efficient swimming strokes through viscous media in two- and three- 
dimensions under various conditions. Marsden (1992) describes a number 
of other examples. 

5 Optics 

Optics has provided fertile ground for studying the geometric phase. (For a 
recent comprehensive survey of experimental activity, see Bhandari (1996)). 
Most investigations have concentrated on polarization effects. The polar- 
ization of a beam of monochromatic, coherent light behaves like a quantum 
mechanical spin, and under cyclic evolutions (and noncyclic ones as well) ac- 
quires a geometric phase similar to that of a spinning particle in a magnetic 
field (cf Section 2.3) There are various ways to change polarization. One 
is by varying the direction of propagation, as the polarization vector must 
remain transverse to it (Tomita and Chaio 1986); in this case the polariza- 
tion has a spin-1, or vectorial character. Another is through optical activity 
induced by anisotropic media, in which case the polarization has a spin- 1/2 
character (Pancharatnam 1956, Bhandari and Samuel 1988, Bhandari 1988). 

A beam of light passing through a cycle of polarization states may be re- 
combined with a second beam whose polarization has remained unchanged. 
Their total phase difference (dynamical + geometrical) can be determined, 
for example, from their interference pattern. In certain cases, the dynamical 
phase shift can be made to vanish. In others, it can be explicitly computed 
and subtracted from the total phase. In either case, one obtains a measure- 
ment of the geometric phase. 

5.1 Geometric phase of coiled light 

In a straight, single-mode optical fibre of circular cross-section, there is 
no preferred polarization; whatever the polarization state of the source, 
it is preserved along the fibre. This is no longer so if the fibre is made 
to bend and wind; then the polarization must change in order to remain 
transverse to the axis of propagation t. Provided the radius of curvature 
R = ||dt/ds|| -1 is large compared to the wavelength A (here s measures 
distance along the fibre), then the polarization vector d(s) of the electric 
displacement D(s,p) = Ree l ( 27rs / A_a; *)/(p)d(s) (where f(p) is the radial 
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modal profile) evolves along the fibre according to 

The corrections to this approximate solution of Maxwell's equations are of 
order X/R. 

For a general elliptical polarization, d is complex. Within the adiabatic 
approximation, both its real and imaginary parts d r and dj remain orthog- 
onal to the direction of propagation t, and each satisfies (5.1) separately. 
As a function of s, i(s) describes a curve C on the unit sphere, and the 
vectors d r (s) and dj(s) are tangent to the sphere at t(s). Equation (5.1) is 
the condition for their parallel transport. (Infinitesimal parallel transport of 
a tangent vector d at t to a tangent vector d + dd at t + dt is accomplished 
by rigidly transporting d to t + dt, and then projecting the result onto the 
tangent plane there.) If the curve C is made to close on itself, the vectors d r 
and dj return rotated through an angle 9 g equal to the solid angle subtended 
by C. 6 g is the geometric phase, depending only on the geometry of C. 

Tomita and Chiao (1986) observed this effect by introducing linearly 
polarized light from a He-Ne laser into a length of helically wound optical 
fibre whose initial and final directions coincided. The direction of linear 
polarization of the emergent beam was measured and found to have turned 
through an angle equal to the solid angle subtended by the cycle of direc- 
tions. The experiment was performed with helices of both fixed and variable 
pitch, as well as for fibres constrained to lie in the xy plane. In this last 
case, the cycle of directions lies on the equator of the unit sphere, and the 
subtended solid angle is a multiple of 2ir. As predicted, the direction of its 
linear polarization remains unchanged. 

The original analysis of the Tomita-Chiao experiment was based on quan- 
tum mechanical considerations, treating a photon propagating in the fibre 
as a spin-one particle whose helicity s.t is constrained to be ±1. (Classically 
this corresponds to left- and right-circularly polarized light.) Just as a spin 
with nonzero magnetic moment tracks the direction of a slowly changing 
magnetic field (cf Section 2.3), the helicity tracks the direction of propaga- 
tion t, and acquires a geometric phase equal to the solid angle described by 
t. By decomposing a linear polarization into a superposition of left- and 
right-circular polarizations, the observed rotation of linear polaration may 
be recovered. 

A systematic derivation, however, requires more detailed analysis, and 
this is most easily carried out within the context of classical (rather than 
quantum) electromagnetic theory. Indeed, the essential features of such an 
analysis are contained in early work by SM Rytov (1938) and VV Vladimirsky 
(1941), who derived and interpreted the parallel transport law (5.1) for po- 
larization in the geometric optics limit (see (Berry 1990) for a discussion). 
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Strictly speaking, geometrical optics is not appropriate for single-mode fi- 
bres (the transverse wavelength is comparable to the fibre diameter). Berry 
(1987a) has carried out a classical analysis using the full Maxwell equations, 
and it is the conclusions of this analysis which are presented above. 

Kitano et al. (1987) performed an experiment related to Tomita and 
Chiao's, in which the cylic change in the direction of propagation was ac- 
complished nonadiabatically through a scries of discrete reflections from 
mirrors. These reflections reverse helicity instead of conserving it, and the 
calculation of the geometric phase must be appropriately modified. 

5.2 Pancharatnam's phase 

When light passes through a homogeneous medium, its direction of propaga- 
tion remains constant. However, if the medium is optically anisotropic, the 
polarization can change. Such changes in polarization are accompanied by 
a geometric phase, discovered in this context by S. Pancharatnam in 1956. 
(See (Berry 1987b) for a contemporary discussion of Pancharatnam's work.) 
For definiteness we consider plane electromagnetic waves propagating in the 
z direction with fixed wavenumber k and frequency oj. The electric field 
E(r,t) = Re y/le l ^ kz ~ ujt ^ \ A) is described by the intensity I a and normalized 
polarization vector | j4) = (A x ,A y ), whose complex components determine 
the phase and relative intensity along x and y. (Dirac notation is used for 
complex two- vectors to highlight an analogy below with two-state systems.) 

Given two such waves, what is their relative phase? Unless their polar- 
ization states | A) and \B) happen to be proportional, it is not clear that the 
question is a sensible one. Pancharatnam found a natural and physically 
relevant prescription for their phase difference by considering the intensity 
Ia+b of their superposition, \/Ta\A) + ^/TB\B). This is given by 

I A+ B = lA + lB + 2^/lA~I^Re(A\B). (5.2) 

For fixed individual intensity I a and Ib, the superposed intensity depends 
only on the interference term, and is maximized when (A\B) = 1. This is 
Pancharatnam's criterion for polarizations being in phase. More generally, 
the relative phase between two polarizations is taken to be arg(^4|i?) (this is 
well-defined unless |^4) and \B) are orthogonal). It is the amount by which 
the phase of one polarization must be retarded relative to the second in 
order that the intensity of their superposition be a maximum. 

Being in phase is not a transitive relation; if \A) is in phase with \B), 
and \B) in phase with |C), then, in general, \C) is not in phase with |j4). 
The nontrivial relative phase axg{C\A) has a simple geometrical interpreta- 
tion. To every polarization state \A) we associate a vector on the unit 
sphere, namely its spin expectation value (^4|cr| A) (a unit vector) rotated 
through ir/2 about the y-axis. This is how polarization is represented, up 
to an overall phase factor, on the Poincare sphere. The y-rotation, a matter 
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of convention, puts the left and right circular polarizations at the north and 
south pole of the Poincare sphere, and the linear polarizations at the equa- 
tor. (In more mathematical language, the Poincare sphere is the projective 
Hilbert space for two-state systems - cf Section 3.3.) 

One can show that the relative phase arg(C|^4) is equal to minus half 
the solid angle subtended by the spherical triangle e^e^ec on the Poincare 
sphere, constructed by joining A, B and C with arcs of great circles. A direct 
demonstration (Pancharatnam's) uses formulas from spherical trigonome- 
try. Another way is to realize that a sequence of polarization states \ip(t)) 
parallel-transported according to (3.14) along great circles on the Poincare 
sphere are, according to Pancharatnam's criterion, in phase with each other. 
Thus arg(C\A) is precisely the geometric phase accumulated under parallel 
transport round e^esec, for which the half-solid-angle formula is discussed 
in Section 2.3. 

To establish a connection between the geometry and an observable effect, 
an analysis of the physical polarization dynamics is required. Pancharatnam 
considered systems of ideal polarizers, or analyzers, which completely trans- 
mit a polarization state and completely absorb the state — orthogonal 
to it in such a way that the incident and transmitted states are always 
in phase. Ideal polarizers are nearly approximated by very thin polariz- 
ing plates. (Mathematically they are represented the hermitian projection 
operator Therefore, if a polarization cycle is generated by ideal 

polarizers, one may conclude that the initial and final polarizations differ by 
Pancharatnam's phase. 

Polarization cycles can also be generated by ideal retarders; these are 
perfectly transparent devices which introduce a relative phase difference be- 
tween two given orthogonal polarization states and — e^. A particular 
example is a quarter-wave plate, for which the orthogonal states are lin- 
early polarized and the phase difference is n/2. (Mathematically, they are 
represented by exp(— ie a ■ cr/2), a two-dimensional unitary matrix with de- 
terminant one.) Polarization cycles generated by ideal retarders produce 
a dynamical phase in addition to Pancharatnam's phase. The dynamical 
phase is easily calculated, however. 

Finally, one can consider optical activity generated by a transparent, 
birefringent gyrotropic medium, which generates relative phase differences 
between orthogonal polarization states whose orientation varies with posi- 
tion. When the variation is sufficiently slow over a wavelength, a beam 
remains in a polarization state determined by the locally uniform dielectric 
tensor. Berry (1987b) applied an adiabatic analysis to Maxwell's equations 
for this system and calculated the phase for a polarization cycle generated 
by a spatially periodic medium. The underlying structure is identical to a 
spin-1/2 particle in a slowly changing magnetic field (Section 2.3). There is 
a dynamical phase and a geometric phase, and the geometric phase is just 
Pancharatnam's. 



44 



Observations of Pancharatnam's phase are based on two-beam interfer- 
ometry. Such experiments were initiated by Pancharatnam himself, though 
current versions are greatly improved with the availability of laser sources. 
Schematically, a plane wave with polarization \A) is split into two coherent 
beams which propagate along the arms of an interferometer. Beam 1 prop- 
agates freely, while beam 2 passes through a sequence of optical elements 
which induce a polarization cycle. The beams are then recombined, and 
their total relative phase determined from their interference. Differences in 
dynamical phases are determined and subtracted, and the geometric phase 
subsequently measured. 

Bhandari (1988) uses laser interferometry methods to measure changes 
in the geometric phase, rather than the geometric phase itself. These are 
readily generated (while keeping the relative dynamical phases fixed) by 
continuously rotating a half-wave plate in the polarization cycle. 

Kwiat and Chiao (1991) observed Pancharatnam's phase in single-photon 
interference experiments using a source of two-photon entangled states. One 
photon, the idler, passed through a narrow-band energy filter to a detec- 
tor. Observation of this photon fixed the energy of its partner, which was 
made to pass through an interferometer as above. Detection rates varied in 
accordance with the intensity, which depends directly on Pancharatnam's 
phase. Various coincidence detectors were used to confirm that genuine 
single-photon interference was being observed. 

Finally, Berry and Klein (1996) perform a modern version of Pancharat- 
nam's experiment, measuring relative phase by direct observation of the 
interference fringes. They also investigate the relation between discrete and 
continuous polarization cycles by varying the number of optical components. 

6 Molecular and Condensed Matter Physics 

Several applications of the geometric phase to molecular spectroscopy (Sec- 
tion 6.1), nuclear magnetic resonance (Section 6.2), and the quantum Hall 
effect (Section 6.3) are discussed. Amongst a number of applications which 
have been omitted from discussion, we mention A-doubling of rotational 
levels in diatomic molecules (Moody et al. 1986), and atomic and reactive 
molecular scattering (Mead and Truhlar 1979, Zygelman 1987, Mead 1992, 
Kuppermann and Wu 1993). 

6.1 Pseudorotation in triatomic molecules 

Ordinarily, the angular degrees of freedom of a spinless system have integer 
quantum numbers. This is because, ordinarily, the angular momentum 
conjugate to an angle 4> is taken to be —ihd^; its eigenfunctions exp(iz^) 
are single-valued only if v is an integer. The Aharonov-Bohm effect (cf 
Section 2.2) provides an exception. In the presence of a magnetic flux line 
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on the z axis, with flux ahc/e, the azimuthal angular momentum of an 
electron is given by —ihd^+a, the constant term a coming from the magnetic 
vector potential A^ s \p, <f>, z) = a/pcj). Therefore, the eigenvalues of are 
no longer integral multiples of h, but are given by instead by (m + a)H, 
where m is integral. 

As discussed in Section 2.7, there emerges in the Born-Oppenheimer ap- 
proximation a molecular analogue of the Aharonov-Bohm effect. For a given 
electric state \n), the nuclear configuration space is threaded by a magnetic 
flux concentrated at points where |n) is degenerate. This flux appears in 
the nuclear Hamiltonian as a geometric vector potential (2.4). As discussed 
by Mead (1992), a triatomic molecule A3 composed of identical atoms (eg, 
alkali metal trimers such as Na3) provides a simple class of examples. A 
striking consequence of the molecular Aharonov-Bohm effect is the appear- 
ance of anomolous half-integer quantum numbers. Their existence has now 
been confirmed experimentally (Delacratez et al. 1986). While this phe- 
nomenon had been antipated by Longuet-Higgens et al. (1958), it is most 
naturally understood in the context of the geometric phase. 

The nuclear configuration space for a triatomic molecule has nine degrees 
of freedom. Six correspond to overall translational and rotational motion, 
which can be ignored in this discussion, leaving three internal vibrational 
degrees of freedom. These may be parameterized by the internuclear dis- 
tances i?i2, -^23) Ri3- Thus, the internal nuclear configuration space Q is 
the space of noncongruent triangles (the i?ij's determining the lengths of 
the triangle's sides) with vertices distinguished. 

In the Born-Oppenheimer approximation, the electronic energy levels 
and eigenstates are regarded as functions on the internal configuration space 
Q. Assuming the molecular Hamiltonian has time-reversal symmetry, we ex- 
pect degeneracies of the electronic states to lie along one-dimensional curves 
in Q (cf Section 2.5). 

It turns out that some of these degeneracies can be determined by sym- 
metry considerations alone. For equilateral configurations R12 = R23 = 
i?i3 = a, the nuclear configuration is invariant under the dihedral group D3 
of symmetries of the equilateral triangle. Therefore, the electronic eigen- 
states \n(a, a, a)) must transform according to an irreducible representation, 
or irrep, of D3 (Tinkham 1964). The group D3 has two one-dimensional 
irreps, associated with states which are either invariant or else change sign 
under the symmetry operations, and a two-dimensional irrep, associated 
with a degenerate pair of states which transform into each other under the 
symmetry operations. Therefore, if \p(a, a, a)) transforms according to the 
two-dimensional irrep, then \p(R\2, R23, R13)) is necessarily degenerate along 
the curve Ru = R23 = R13 = a (with a > 0) in Q. Under parallel-transport 
round this curve (cf (2.2)), \p) changes sign. 

One might expect the energy E P {R\2, R23, R13) of this state to assume 
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its minimum in an equilateral configuration. The Jahn- Teller effect (Jahn 
and Teller 1937) shows this is not the case. It turns out to be energetically 
favourable for the molecule to assume an asymmetrical configuration, so 
that when it is in its p-electronic state, the ground vibrational state of the 
molecule is distorted. 

To analyze the vibrational motion, it is useful to introduce cylindri- 
cal coordinates (p,(p,z) in Q. The vertical coordinate z, given by R\ 2 + 
i?23 + #13, determines the length scale; oscillations in z describe uniform 
dilations and contractions of the nuclei. The radial coordinate p, given 
by 2R 2 12 + (V3 - - (\/3 + l)i?f 3 , describes the de gree of distortion; 

p = corresponds to equilateral configurations. Thus, it is along the z axis 
that the electronic state \p(p, (j), z)) is degenerate. The azimuthal angle (f), 
given by tan _1 [v / 3(i?23 ~~ ^13)/ (^12 ~ #23 ~~ #13)]) describes the direction 
of the distortion. Displacements in (f> are called pseudorotations. A 2ir- 
pseudorotation returns the nuclei to their original positions, while a 2ir/2> 
rotation induces a cyclic permutation amongst them. 

A systematic calculation of the vibrational Hamiltonian H vi ^ requires 
an explicit expression for the electronic energy E p (R±2, R23, R13) along with 
a careful derivation of the kinetic energy in the (p, <fi, z) coordinates. Such 
a treatment is described by Mead (1992); however, for this discussion an 
approximate form will suffice. The potential is taken to be harmonic in z 
and p and independent of (j), and the kinetic energy to be of the standard 
form. Thus 

= m + 4? + ii + 5 M "> - "» )2 + l M ^ z - z " )2 ' (<u) 

where (po,(f>o) denotes the asymmetrical configuration which minimizes E p . 
The vibrational energy levels are characterized by three quantum numbers 
j, m, k, for the respecive degrees of freedom p, (f>, z. They are given approxi- 
mately by 

1 (mh) 2 1 

E(j, m, k) = (j + -)hu p + ^f- + {k + -)hw z , (6.2) 

where the azimuthal moment of inertia / is taken to be Mp^ (the approxima- 
tion consists of ignoring the dependence of / on the radial quantum number 
j)- 

The presence of the electronic degeneracy along the z axis introduces 
a geometric vector potential A p = a<p/p into the nuclear Hamiltonian, 
which accounts for the sign change of the electronic eigenstate under parallel 
transport. (In the chosen coordinates, a = |.) Its effect is to replace the 
pseudoangular momentum operator = —ihd^ by —ihd^ + just as in 
the magnetic Aharonov-Bohm effect. Thus, the pseudorotational quantum 
number m in (6.2) is half-integral. 
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The half-integral quantization of pseudorotation was confirmed in a sys- 
tematic spectral analysis of the vibrational spectrum of Na3 be G Delacretez, 
ER Grant, RL Whetten, L Woste, and JW Zwanziger (1986). A sequence 
of rovibrational levels corresponding to radial and pseudorotational motion 
were identified. It was found that the observed spectrum could be fitted to 
(6.2) only by assigning half-integer quantum numbers to the pseudoangu- 
lar momentum. The result is an excellent agreement between observed and 
calculated energy levels. 

6.2 Nuclear magnetic resonance 

An rf magnetic field applied to polarized spins generates transitions between 
their stationary states. When the field is in resonance with the spin preces- 
sion frequency, the transition probabilities can become large, even though 
the rf field may be weak. Transitions are manifested in oscillations of the spin 
magnetic moment, whose amplitude and frequency can be measured with 
great accuracy. This is the basis for nuclear spin resonance experiments. 
In nuclear magnetic resonance (NMR), the nuclear spins are polarized by 
a static magnetic field. In nuclear quadrupole resonance (NQR), the spins 
are polarized by an electric field gradient (in the ambient crytalline poten- 
tial, for example) which couples to the spin though the nuclear quadrupole 
moment. 

Moody et al. (1986) and Cina (1986) pointed out that cyclic variations 
in the rf field would generate geometric phases in the nuclear spins, which 
could be detected as shifts in the frequency spectrum of the nuclear magnetic 
moment (this is an example of the third type of experiment described in 
Section 2.6). Pursuing this suggestion, several groups have carried out NMR 
studies of the geometric phase, including its adiabatic (Tycho 1987, Suter 
et al. 1987), nonadiabatic (Suter et al. 1988), and nonabelian (Zwanziger et 
al. 1990b) versions. These experiments have not only confirmed theoretical 
predictions, but have also brought to light interesting new phenomena in 
NMR spectroscopy. 

In an NMR experiment, the nuclear spins (let s = 1/2, say) are initially 
polarized by a strong static magnetic field B z z; the spin eigenstates |±) 
have energies ±hwo/2, where oj$ = jB z is the Larmor frequency and 7 
the gyromagnetic ratio. For nonstationary states, the magnetic moment 
M(t) = h / y(ip(t)\a\ip(t)) /2 precesses about z at the Larmor frequency, so 
that its frequency spectrum has a single peak at ujq. 

A weak, circularly polarized rf field, B r f(i) = B x cos wtx + BySinuty, 
is then applied. (Weak means that _B r f <C B z ). The subsequent dynamics 
is most easily analyzed in a rotating frame turning about z with angu- 
lar frequency oj. There the effective magnetic field B rot is static, and the 
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Hamiltonian is given by 

#rot = ^7^ B rot • cr, B rot = {B x ± + B y y + (B z - u/j)z) . (6.3) 

The spin eigenstates |±) rot in the rotating frame are polarized along B rot 
with energies ±foj rot /2, and the magnetic moment M rot (i) precesses about 
B rot with frequency ui TOt , where u TO t = 7B rot . Near resonance, where uj ps 
ujq, the axis of precession is nearly perpendicular to z (even though B T f <C 
B z ). 

The magnetic moment M(i) in the laboratory frame is obtained by turn- 
ing M rot (£) through an angle — ujt about z. The effect of the rf field is to 
superpose a slow longitudinal nutation (of frequency uj r ) on the fast Larmor 
precession of M(i) about z (of frequency ujq). In consequence, the frequency 
spectrum of the magnetic moment has additional peaks at ujq ± uj TO t- 

The Hamiltonian H ro t of (6.3) depends parametrically on the compo- 
nents of the rf and static magnetic fields, B x , B y and B z , and vanishes for 
B x = By = and B z = w/7. If these components are slowly modulated 
through a cycle C, with a period T much greater than the precession period 
2ir/uj TOt in the rotating frame, then the eigenstates |±) rot of H rot evolve adi- 
abatically, and acquire with each modulation cycle a geometric phase 7±(C I ) 
given by half the solid angle, a, described by C with respect to (0, 0, u/j) (cf 
Section 2.3). These geometric phases produce shifts ^^ha/T in the energies 
±^Huj Iot of the unmodulated Hamiltonian H Iot , and in turn shift the preces- 
sion frequency of the magnetic moment from uj TOt to uj TOt — a/T. (This shift 
is the magnetic analogue of the daily precession of the Foucault pendulum 
discussed in Section 4.1.3.) 

These frequency shifts were observed by D Suter, GC Chingas, RA Har- 
ris, and A Pines (1987) in the NMR response of proton spins in a wa- 
ter/acetone sample. The adiabatic cycle C was produced by placing the 
magnetization detector in a frame rotating with angular frequency oj^ct 
about z. In the detector frame, the effective magnetic field turns about 
z with angular frequency 5 = uj — wjct! f° r small enough S, the spins adia- 
batically track the turning field. A variation of this experiment was subse- 
quently performed by D Suter, KT Mueller, and A Pines (1988), in which 
the geometric phase was observed as a rotation in the spin-echoed magnetic 
moment. This technique enabled observation of the individual geometric 
phases j± rather than the difference 7+ — 7_ . The spins were taken through 
nonadiabatic cycles, thereby incorporating the Aharonov-Anandan phase (cf 
Section 3.3). 

In a nuclear quadrupole resonance experiment, the spin polarization is 
induced by an inhomogeneous electric field rather than a static magnetic 
field. Gradients in the electric field couple to spin through the nuclear 
electric quadrupole moment. The effective spin Hamiltonian is of the form 
Hq = Hujq(S • n) 2 , where S denotes the dimensionless spin operators and n 
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the symmetry axis of the electric field gradient tensor (assumed to be axially 
symmetric). The eigenstates |m(n)) are polarized along n with energies 
hujQim? . Note that these are doubly degenerate for 

Geometric phase effects in the NQR spectrum of 35 CI (for which s = 3/2) 
were studied by R Tycko (1987). In the experiment, a coherent superposition 
of eigenstates of Hq was created at t = by an rf pulse; the induced 
magnetic moment M(i) was then observed and spectrally analyzed. The 
spectrum of Hq consists of two degenerate doublets with energies \Tvjjq 
and \fvjJQ (corresponding to m = ±3/2 and m = ±1/2 respectively), so 
that the magnetic moment spectrum has a single peak at frequency 2u>q. 
Adiabatic cycling of Hq was achieved by rotating a crystalline sample of 
NaC103 with angular frequency ur <C ujq about an axis inclined to the 
symmetry direction n. With each rotation, the spin eigenstates acquire 
geometric phases (3 m (these may be computed explicitly), which shift their 
energies by h(3 m u}R/2ir and lift the double degeneracies. For the particular 
geometry realized in Tycko's experiment, theory predicts that the single 
peak in the magnetic moment spectrum is split into three at frequencies 2uq 
and 2coq ± \/3ujr/it. Precisely this splitting was observed in the experiment. 

Since Hq has degenerate energy levels, one might expect nonabelian 
geometric phases (cf Section 3.2) to appear. As shown by Tycko (1987) and 
subsequently by A Zee (1988), the experiment can be analyzed completely 
in terms of the ordinary (abelian) geometric phase. The nonabelian vector 
potentials A( ±3 /2) and A( ±1 / 2 ) of (3.5) (there is one for each degenerate 
doublet of Hq) can be evaluated explicitly; they involve matrix elements 
of the spin angular momentum between degenerate states. For m = ±3/2, 
the off-diagonal matrix elements vanish identically (S has nonzero matrix 
elements only between states with m differing by one), and the diagonal 
matrix elements yield the ordinary geometric phases P±z/2- For m = ±1/2, 
these off-diagonal elements do not vanish. However, for the adiabatic cycles 
used in the experiment, in which n is turned about a fixed basis for 

the m = ±1/2 doublet can found which diagonalizes A( ±1 / 2 )- During the 
cycle there are no transitions between these basis vectors, and each acquires 
an ordinary geometric phase (3±i/2- 

Zee (1988) shows that more general adiabatic cycles of n generate gen- 
uine nonabelian geometric phases in the NQR spectrum. JW Zwanziger, M 
Koenig and A Pines (1990b) performed an experiment to observe these non- 
abelian effects. The more complicated adiabatic cycles were generated by 
mounting the NaC103 crystal on a double rotor. Careful analysis predicts 
that the NQR resonance should be split into five peaks instead of three. 
This is precisely what was observed, with excellent quantitative agreement 
between theory and experiment. 
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6.3 The quantum Hall effect 

Charged particles in crossed uniform electric and magnetic fields experience 
a net drift along the direction E x B normal to both - this is the basis for the 
Hall effect. The Hall effect concerns specifically the behaviour of electrons 
in a planar conductor to which a perpendicular magnetic field B = Bz and 
a tangential electric field E = Eit are applied. The resulting current density 
has a transverse component j y = ohE. an is the Hall conductance. Using 
classical theory, one obtains for large magnetic fields that 

c nec 

where n is the two-dimensional electron density. 

(Alternatively, one can consider conductors of finite transverse extent, 
eg long thin wires. Then the current density is necessarily longitudinal, 
whereas the electric field has a transverse component. In this case it is easier 
to compute the Hall resistance pu = E y /j. Classical theory gives p c H = 
B/nec for all magnetic fields, not just large ones. The Hall conductance and 
resistance are the off-diagonal elements of the conductance and resistivity 
tensors a and p = a" 1 , and the relation between them, an = Ph/(Po + Ph)i 
depends on the longitudinal or Ohmic resistivity po-) 

For large magnetic fields and at low temperatures, the classical predic- 
tion (6.4) fails; what is actually observed is a spectacular phenomenon, the 
quantum Hall effect. As B is increased, instead of falling off monotoni- 
cally as 1/B, the Hall conductance goes through a series of plateaus where 
its value is constant, separated by intervals of monotonic decrease. These 
constant values are sample- independent, and satisfy with extremely high 
accuracy (about one part in ten million) the quantization condition 

e 2 

a H = r—. (6.5) 
he 

As B increases, the coefficient r takes on integer values, and decreases with 
B in unit steps until it reaches one. This is the regime of the integer quantum 
Hall effect (von Klitzing et al. 1980). As B is increased still further, one 
enters the regime of the fractional quantum Hall effect (Tsui et al. 1982); r 
takes on rational values (< 1) with odd denominators. 

The theory of the quantum Hall effect begins with the quantum me- 
chanics of two-dimensional electrons in magnetic fields. The one-electron 
stationary states are the Landau levels Em = hui c (N + 1/2), N = 0, 1, 2, . . ., 
where ui c = eB/mc is the classical cyclotron frequency. The Landau levels 
are highly degenerate; each can accomodate B/<&o electrons per unit area 
($o = hc/e is the quantum of magnetic flux). In real materials, the crys- 
talline potential, impurities, and electron-electron interactions all serve to 
break this degeneracy. However, for sufficiently large magnetic fields, there 
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(6.4) 



is negligible mixing between the Landau levels, so that each is split by the 
perturbations into B/^q sublevels per unit area. 

If the electron density in the sample is n, then the ground state has the 
first v = n&o/B Landau levels filled, the uppermost only partially if v has 
a fractional part. (Thus, the larger the -B-field, the fewer occupied Landau 
levels.) The integer quantum Hall effect can now be explained by making 
two assumptions: i) only a small proportion of sublevels in each Landau 
level can carry current, and ii) each filled Landau level contributes precisely 
e 2 /he to the Hall conductance. Together, these imply that a nearly empty 
Landau level contributes nothing to the Hall conductance, as the occupied 
sublevels are unlikely to be the current-carrying ones, whereas, conversely, 
a nearly full one contributes e 2 /he. From these qualitative considerations 
follows &h = Me 2 /he {[v\ denotes the integer nearest v), precisely (6.5), 
with r = \y\. 

The first assumption can be justified on the basis of localization the- 
ory - most of the sublevels are localized, whereas only extended states can 
carry current. For the second assumption, DJ Thouless, M Kohmoto, MP 
Nightingale and M den Nijs (TKN 2 ) found, in 1982, a topological expla- 
nation. In retrospect, their argument can be seen as an application of the 
theory of the geometric phase. (Indeed, this connection was drawn as early 
as 1983 by Simon.) 

In outline, their analysis is as follows. We consider an electron in two 
dimensions in a uniform magnetic field Bz and periodic crystalline poten- 
tial V(x,y), with periods a and b in x and y. Because the vector po- 
tential A = (By/2,—Bx/2) is not periodic, neither is the Hamiltonian 
mv 2 /2 + V(x,y), in which v = (p + eA/c)/m is the velocity. However, the 
Hamiltonian is invariant under the magnetic translation operators T x and 
T y , which shift x and y by a and b while leaving the velocity (rather than 
the momentum) unchanged. (Explicitly, T x = ex.p(—iap x /h + iirBay/§o) 
and T y = exp(—ibp y /h — mBxb/^o).) 

While the magnetic translation operators commute with the Hamilto- 
nian, they do not necessarily commute with each other. In general, T x T y = 
exp(27riBab/<&o)T y T x . But if we take the flux Bab through a unit cell to 
be a rational multiple p/q of the magnetic flux quantum $o> then (T x ) q - 
a translation by qa in x - and T y do commute, and Bloch's theorem ap- 
plies. In this case, the energy levels form a band structure, whose bands 
e n (k) are periodic over the first Brillouin zone [0,27r/ga] x [0, 2ir/b\. The 
stationary states V'k,n( r ) are Bloch waves exp(ik.r)«k jn (r), where u^ n (r) is 
an eigenfunction of the Hamiltonian 

H(k) = ^m(-ihV + eA/c + hkf + V(x, y) (6.6) 

with energy ek, n satisfying the generalized periodic boundary conditions 
Uk, n {x,y) = exp(iTrBqay/<& )u^ n (x+qa,y) = exv{-inBxb/§o)u^ n (x 1 y+b) 
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over the unit cell. (These are obeyed rather than ordinary periodic bound- 
ary conditions because we are dealing with magnetic rather than ordinary 
translation operators.) For convenience we will sometimes denote the indices 
re, k by a single label a. 

The Hall conductance can now be calculated perturbatively. An applied 
electric field Ex generates perturbed eigenstates \u' a ) = \u a ) + \Su a ), and 
induces transverse current densities j y = (ThE. To lowest order in pertur- 
bation theory, these are given by — e • 2Re (u a \v y \5 a ) /Aq. Here Aq is the 
area of the sample, and the transverse velocity v y = —ihd y — eBx/c + hk y 
includes a contribution from the Bloch vector. Note that v y may be ex- 
pressed as dH/dhk y . From first-order perturbation theory, (u/3\5u a ) = 
{up\ — eEx\u a ) / (e(a) — e(/3)),/3 / a. Combining these results with the 
identity [x,H] = —ihv x , we obtain the Kubo formula 

Im fe^P (6 ' 7) 

for the Hall conductance of the ground state. It includes contributions from 
all states below the Fermi energy ep. 

Suppose we regard HQs.) as a family of Hamiltonians parameterized by 
the Bloch vector k. Under parallel transport round a closed cycle C in the 
first Brillouin zone, the eigenstate |ti n ,k) acquires a geometric phase 7 n ,k(C)> 
which may be expressed (cf Section 2.2) as the integral of the flux field 
V n (k) = Im (Vk^ n ,k| x I Vk^ n ,k) over the interior of C. Amongst the various 
expressions for V n (k) given in Section 2.2, we note that (2.10), which involves 
matrix elements of the gradient of the Hamiltonian, is contained precisely in 
the formula (6.7) for the Hall conductance. This allows us to rewrite (6.7) 
as oh = e VA)frc]C n !Ck ^ri(k). Finally, for samples large compared with 
the unit cell, the sum over Bloch vectors may be replaced by an integral 
As /Ait 2 J d 2 k over the first Brillouin zone. Thus we obtain 

n 

where a n = J (PkYnk is the conductance (in units of e 2 /he) of the nth band. 

TKN 2 show that a n is an integer. This follows by first noting that V n (k) 
is periodic over the Brillouin zone (even though |n„ 5 k) is not), so that the 
domain of the k-space integral is a closed torus T^. Indeed, one can show that 
Ki(k) is given by Im (VkVVi,k| x | VkVVi,k) phis a perfect curl (whose integral 
over the torus vanishes by Stokes' theorem). The Bloch states IV'n.k) ar e 
periodic (up to a phase factor!) over the first Brillouin zone, and constitute 
a complex line bundle (cf Section 3.1.4) over with curvature V n (k). The 
integral over Tk of the curvature over the torus is the first Chern number, an 
integer. For large magnetic fields, the bands coalesce into nearly degenerate 
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sublevels (p per band) associated with the Landau levels of a pure magnetic 
Hamiltonian. One can show that sum of Chern numbers from a given Landau 
level is equal to one. This completes the topological explanation for the 
quantization of Hall conductance. TKN 2 also calculate the Chern numbers 
for individual sublevels, which are obtained from the solutions of certain 
Diaphontine equations in p and q. 

TKN 2 was the first of a number of important applications of geometric 
phase methods to the study of the quantum Hall effect. Several subsequent 
works address the fractional effect, where electron-electron interactions come 
into play. Niu, Thouless and Wu (1984) show how the results of TKN 2 could 
be generalized to accomodate them. Arovas et al. (1984) calculate properties 
of the elementary excitations of Laughlin's many-body fractional-quantum- 
Hall ground state (Laughlin 1983) using geometric phase arguments. Con- 
cerning the integer case, Arovas et al (1988) establish a topological criterion 
for localization in terms of Chern numbers. Wilkinson (1984) incorporates 
the geometric phase into semiclassical calculations of the single-particle en- 
ergy spectrum. 

Finally, we mention studies by Avron and coworkers (Avron et al. 1988, 
Avron 1995) of electron transport in quantum systems. While not specifi- 
cally concerned with the quantum Hall effect, they elucidate essentially sim- 
ilar topological mechanisms. One class of models consists of planar networks 
of one-dimensional wires whose closed loops Cj are threaded by magnetic 
fluxes 4>i normal to the plane. The physical properties of the network are <3?o- 
periodic in the fluxes (fluxes differing by an integral number of flux quanta 
are related by a gauge transformation). Fixing all but two of the fluxes, 
the stationary state \ip n ((pi,(pj)) constitute a complex line bundle over the 
two-torus of fluxes (the analogue of the first Brillouin zone) . The associated 
Chern number, the integral of the curvature, has the following physical in- 
terpretation: it represents the average with respect to fa of the number of 
electrons transported round Q as <pj is increased adiabatically by one flux 
quantum. 

Glossary 

dynamical phase The time integral of the energy of an evolving quantum 
state divided by Planck's constant h. 

geometric phase A phase factor accumulated by a cycled quantum eigen- 
state which depends only on the geometry of its evolution. 

geometric magnetism Effective Lorentz-like force which appears in the 
dynamics of a slow system coupled to a fast one after the fast degrees 
of freedom have been averaged away. 
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Hannay angle Displacement in the angle variables of adiabatically cycled 
integrablesystems of purely geometric origin. 

holonomy Change in phase induced by parallel transport of a quantum 
state around a closed curve. 

integrable system A Hamiltonian system with as many constants of the 
motion in involution as there are degrees of freedom. 

parallel transport A prescription to determine the overall phase of an 
evolving quantum state. 

reduction General theory of Hamiltonian dynamics in the presence of sym- 
metry. 
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